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Abstract: Several observables for the deeply virtual Compton scattering process have 
4^ been simulated in the kinematic regime of a proposed Electron-Ion Collider to explore 

the possible impact of such measurements for the phenomenological access of generalized 

^ parton distributions. In particular, emphasis is given to the transverse distribution of sea 

^ quarks and gluons and how such measurements can provide information on the angular 

momentum sum rule. The exact lepton energy loss dependence for the unpolarized t- 
— jl differential electroproduction cross section, needed for a Rosenbluth separation, is also 

reported. 
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1 Introduction 

During the last decade the collaborations at the Hadron Electron Ring Accelerator (HERA) 
and the Thomas Jefferson National Accelerator Facility (JLAB) spent lately significant 
effort to measure exclusive processes such as the electroproduction of a real photon (a 
process known as deeply virtual Compton scattering (DVCS) [1-18]), vector mesons (VM) 
p° [19-29], (j) [24, 30-35], uo [36, 37], J/ip [20, 38-40], T [41-43], and the pseudoscalar meson 
7r + [44-46] in the deeply virtual region in which the virtuality Q 2 > 1 GeV 2 of the exchanged 
space-like photon allows to resolve the internal structure of the proton. The HERA collider 
experiments [1-6]. found that the exclusive cross sections grow with increasing energy W, 
where the effective "pomeron" intercept is larger and the slope parameter smaller than for 
the soft pomeron trajectory [47], introduced to describe elastic (anti-)proton- proton high 
energy scattering. Moreover, the exponential t-slope parameter as a function of the scale 
Q 2 + My M was determined by fitting the ^-dependence of the cross section for exclusive 
vector meson production and DVCS [43], which makes loose contact to the idea of imaging 
the proton content [48]. 

Various phenomenological and theoretical descriptions for these exclusive processes 
have been proposed and utilized. In the high-energy region it is popular to understand these 
processes in terms of the pomeron picture [49], perturbative high-energy QCD [50, 51], the 
color dipole picture [52, 53], or in terms of the color glass condensate approach [54, 55]. 
In the deeply virtual regime exclusive processes provide an important tool in accessing the 
generalized parton distributions (GPDs) [56-58], bridging thereby the high and medium 
energy regions. GPDs also enter in the hand bag model approach [59-61], which allows 
to describe observables that in the perturbative GPD approach are considered as non- 
factorizable contributions that cannot be perturbatively treated. In all these approaches 
the underlying mechanism is a ^-channel exchange with different degrees of freedom. 

Based on factorization theorems [62, 63], GPDs offer a partonic interpretation of these 
processes, where unobserved transverse degrees of freedom are integrated out. Thereby, 
these universal functions, defined in terms of matrix elements of quark and gluon operators 
or, alternatively, as a non-diagonal overlap of light-cone wave functions [60, 64, 65], encode 
the non-perturbative aspects of the nucleon. Because of their fundamental QCD definition, 
a whole framework is built up around GPDs, various aspects of the GPD framework are 
reviewed in [66, 67]. In particular, GPDs provide an access to the transverse spatial 
distribution of patrons [68-70] , and appear in the gauge invariant decomposition of the 
nucleon spin in terms of quark and gluon degrees of freedom [71]. 

Phenomenologically, exclusive electroproduction of a real photon, DVCS, diagrammat- 
ically depicted in Fig. 1 (left), is the golden channel to constrain GPDs as it is theoret- 
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Figure 1. Amplitudes contributing to the photon leptoproduction cross section in leading order 
approximation of QED: the DVCS amplitude (left) while the remaining two diagrams (middle and 
right) represent the Bethe-Heitler amplitudes, parameterized by hadronic electromagnetic form 
factors. 

ically clean and the phase of its amplitude can be measured using the interference with 
the Bethe-Heitler (BH) amplitude (see Fig. 1 middle/right). Besides that, the measure- 
ment of Compton scattering observables, even at rather low photon virtuality, is important 
since it provides insight into the fundamental Compton scattering process in the virtual 
regime. Since the virtual Compton process contains twelve helicity amplitudes (or equiva- 
lently twelve complex Compton form factors (CFFs) [72]), their disentanglement is already 
an experimental challenge. The measurement of CFFs should be considered a primary 
task, as important as the measurement of electromagnetic nucleon form factors. In return, 
the (partial) disentanglement of the various CFFs offers then a phenomenologically much 
simpler and cleaner access to GPDs. 

Based on present phenomenological GPD knowledge, Monte Carlo simulations, and 
GPD fitting routines, we explore in our studies here both the DVCS process and the 
access to the spatial transverse distribution of quarks and gluons at a proposed Electron- 
Ion Collider (EIC). The much more general physics case of this suggested high-luminosity 
collider with a dedicated detector for exclusive channels in the medium to high energy 
regime of lepton- nucleon and lepton- nuclei scattering is described in [73] . 

The rest of this article is organized as follows: in Sect. 2 we introduce the theory 
elements, needed for the access of GPDs from DVCS observables, including also, for the 
unpolarized case, the exact dependence of a (reduced) ^-differential photon electroproduc- 
tion cross section on the electron energy loss variable y. Furthermore, we give a short 
overview of existing DVCS measurements. In Sect. 3 we describe the planned EIC at its 
different stages and the Monte Carlo simulation technique used in the generation of EIC 



DVCS pseudo-data. In Sect. 4 we shortly introduce three GPD models which are then 
utilized to provide predictions for the ^-differential DVCS cross section, single spin and 
lepton charge DVCS asymmetries at different EIC kinematics. In Sect. 5 we discuss the 
access of GPD H and E at the final stage of EIC by using the DVCS cross section and 
single transverse proton spin asymmetry. Furthermore, we quantify the implications of 
such measurements for the imaging of the proton and comment on the qualitative aspects 
of such measurements for the spin sum rule. Finally, we summarize and conclude in Sect. 6. 

2 Deeply virtual Compton scattering 

The differential photon electroproduction cross section is five-fold and consists of the sum 
of the BH amplitude squared, DVCS amplitude squared, and the interference (INT) terms, 
where the latter is charge odd: 

dcr ep^e P1 ^ da e ^ e ^' BH (Fi, F 2 ) da e ^ e ^' INT (Fi, F 2 , J 7 ) d cr ep^ep7,(D)vcs^ 5 jr*) 



dxB did Q 2 d(j)d(f dx^dtdQ 2 d(j)d(f dx^dtdQ 2 dc/>dip dx^dtdQ 2 d(j)dLp 

(2.1) 

Here the +(— ) sign is valid for electron (positron) beam, xb is the common Bjorken scaling 

variable, <j) is the azimuthal angle between lepton and hadron scattering planes, and (p = 

$ — 0, where <£>(= (/)$) is the angle between the lepton scattering plane and a possible 

transverse spin component of the incoming proton at rest. We adopt in the following 

the frame conventions of [72] (virtual photon momentum is counter-along the z-direction 

and x-component of the incoming electron momentum is positive). To the leading order 

(LO) in the electromagnetic fine structure constant a em = |^ « t-U, and neglecting the 

electron mass, the three terms on the r.h.s. of (2.1) are exactly known in terms of the 

electromagnetic Pauli form factor F\(t) and the Dirac form factor i 7 ^), parameterizing 

the BH amplitude, and a set of twelve photon helicity dependent CFFs ^(xb,^ Q 2 ), 

parameterizing the DVCS amplitude, see Fig. 1. These CFFs are labeled by the helicities 

of the incoming a E {+, 0, — } and outgoing photon b E {+, — } and they are called 

Tab £ {Hah, £ab, Hab, £ab} with Jq- = ^0+ j ^+- = ^-+ ? (2.2) 

more details can be found in [74, 75]. Analogously to the Dirac and Pauli form factor F\ 
and F2 (axial and pseudo-scalar form factors Fa and Fp), the CFFs %{%) and £{£) ) are 
associated with conserved proton helicity amplitudes and helicity flipped ones, respectively. 

The three separate terms of the differential electroproduction cross section (2.1) can 
be expanded w.r.t. harmonics of the azimuthal angle 0, 

da e P ^n {FuF2) _ aL , B l(l + e 2 )- 5 / 2 f * 1 

,n=0 ) 

(2.3) 



dx B dtdQ 2 d<pd<p 16tt 2 Q 4 tVi{4>,y)V 2 {(t>,y) 



j ep^em,mT/p p\ 3 -lfi , 2\-l/2 ( 3 ^| 

^- ,^ 9 \ I = ,. 7U, I/,'!, , U NT + y [c™ T cos(n</>) + 4 NT sm(n<f>)} I , 

(2.4) 

(2.5) 

Here l/('Pi(<f>,y)'P2{ ( t>,y)) are (rescaled) BH propagators, defined in (32) of [72], the energy 

ioss 

1 Q 2 

V= -^J=M$ (2 - 6) 

of the electron depends for fixed xb and Q 2 on the center-of-mass (c.o.m.) energy squared s, 
and, finally, we used the shorthand e = 2x^M p /Q. Moreover, all of the Fourier coefficients 
c^ and s'n depend on the polarization vectors of the protons. The explicit expressions for 
an incoming polarized nucleon has been presented in [75] , where for transverse polarization 
the coefficients can be further decomposed in cos((/?) and sin((/?) harmonics. Note that if 
one reduces the five- fold cross section (2.1) to a four- fold one by integrating over (/?, the 
(/^-harmonics drop out and the remaining unpolarized and longitudinally polarized parts of 
the expressions (2.3-2.5) are multiplied by a factor 2tt. The knowledge of the coefficients 
in the BH term (2.3) is limited only by the knowledge of the proton form factors F\(i) and 
F2(t). The coefficients of the interference (2.4) and (D)VCS term (2.5) are linear and bi- 
linear in the CFFs, respectively. We emphasize that electromagnetic corrections will enter 
in all three terms. So far such a em /^-proportional corrections are only partially taken into 
account in radiative correction procedures. 

Adopting the discussion of [72], we can state that an over-complete set of observables 
exist and that at least in principle their experimental measurements would allow to extract 
the real and imaginary parts of all twelve CFFs (2.2). Loosely speaking, in the deeply 
virtual regime the first harmonics in the interference term are dominant, i.e., proportional 
to 1/Q 3 , and are governed by twist-two associated CFFs (or GPDs), while the constant 
and second harmonics are kinematical suppressed by 1/Q and arise in leading order of 
perturbative QCD from both twist-two and twist-three associated CFFs (or GPDs). The 
third harmonics are counted as leading twist contributions, however, they arise in next- 
to-leading order (NLO) of perturbative QCD from gluon transversity GPDs. A rather 
analogous counting scheme holds for the zeroth, first, and second harmonics of the DVCS 
term, where, compared to the interference term, an additional kinematical factor 1/Q 
appears. Because of this mismatch in twist and power counting, some care is needed. 

In the rest of this Sect. 2 we consider: in Sect. 2.1, the ^/-dependence of the ^-integrated 
electroproduction cross section (2.1) for an unpolarized proton and in Sect. 2.2 we point out 



that the relation of helicity CFFs to GPDs can be systematically improved. In Sect. 2.3 
we give a short overview of existing DVCS measurements and make a loose contact to 
CFF/GPD phenomenology. 

2.1 Rosenbluth separation of electroproduction cross section 

It would be very desirable to decompose the photon electroproduction cross section (2.1) 
into its different parts (2.3-2.5). In an experimental setup in which both electrons and 
positrons are available, the charge-odd interference term (2.4) and the charge-even part, 
given as sum of BH and DVCS cross sections (2.3,2.5), can be obviously separated from each 
other by forming the difference and sum of electron and positron cross sections. Having 
only an electron beam at hand, it remains so far unclear to what extent a variation of 
c.o.m. energy (or electron/proton beam energy) allows for a Rosenbluth separation, which 
is expected to be much more intricate than in the case of elastic form factors or deeply 
inelastic scattering (DIS) structure functions. We recall that in these cases two form factor 
combinations (or structure functions) enter the unpolarized cross sections; however, both 
of them arise from transversely or longitudinally polarized photon exchanges and are thus 
accompanied with a different y(xB, Q 2 ,s) dependence, which varies for fixed xb and Q 2 
with the c.o.m. energy y/s, see (2.6), i.e., with the beam energy (or energies). 

Having the exact analytic expressions of [75] in mind, it looks hopeless to employ a 
Rosenbluth separation directly to the five- fold (or four- fold) cross section (2.1). Thus, it 
is more appropriate to project first on the azimuthal angle harmonics, where, however, 
the <j)- and ^/-dependencies of the BH propagators should be treated in such a way that 
the final result is most appropriate for the analyzes of experimental data. Including these 
propagators in the integral, as done in [72], provides a truncated Fourier series and allows 
for a rather simple power counting scheme; however, these Fourier coefficients will not have 
a simple y-dependence. Alternatively, one may stick to the standard Fourier coefficients, 
e.g., calculated from 

#cos(n< ^ B ^ Q2# for „ = 0,1,2,3,.- , (2.7) 

where the DVCS cross section only contributes to the first three lowest coefficients. In 
the following we will first consider only the lowest harmonic, i.e., n = 0, which has a 
surprisingly simple and obvious ^/-dependence. 

Let us first introduce a formally defined ^-differential "photoproduction" cross section. 
It is obtained by integrating the four-fold electroproduction cross section over the azimuthal 
angle 0, multiplying it with an infinitesimal electron phase space element, and dividing it 



by a flux factor, 

da TOT (x B ,t,Q 2 \y) _ i r ^d* ep ^ ep i(x B ,t,Q 2 \y) , _ 2 



dt r(xB, Q 2 |?/) 7-7T ' dtd(f)dxBdQ 2 



/IT 
-7T 



<fe B dQ (2.6 



= dc7 BH (x B ,t,Q 2 |i/) ± dc7 INT (xB,t,Q 2 b) | dc7 DVCS (x B ,t,Q 2 |y) 

dt dt dt 

where as before the positive (negative) sign of the interference term refers to an electron 
(positron) beam. For the virtual photon flux we adopt the Hand convention [76] by taking 

r(x B ,Q 2 b) = — y \^~*o with e(y) = 1 ~ y ~ L ^~ 22 (2.9) 

where e(y) is the ratio of longitudinal and transverse photon flux. 

The ^/-dependence of the three terms in (2.8) and the explicit expressions for the 
Fourier coefficients can be evaluated from (2.3-2.5) . Thereby, the DVCS cross section 
is the most simplest one and given by the constant harmonic in (2.5), which is further 
specified in (36) of [75]. The BH cross section can also be analytically calculated, where 
due to the ^-dependence of the BH propagators also higher (f) harmonics that arise from the 
interference of photon helicity flip amplitudes, specified in (35-37) of [72], enter. Thereby, 
the integration over the azimuthal angle <j) generates a characteristic y-dependent function 
that stems from the product l/{Vi((j),y)V2{4>,y)) of BH propagators. Consequently, this 
function inherits the i^-channel pole of one BH propagator at 

y = y C o\ with y co i = — 2- — - , 

Q z + x B t 

where the real photon and incoming electron momenta are collinear. In the following we 
present results for the region y < y co i, in which this characteristic function reads: 

/\ | XBt 

1 _ V Q 2 



Q 2 y i - . Q 



1 + Si - vM l ~ x b) ! + ^-?/ (! + !#) 1 + Ss ~ vM l ~ x b) [y C oi - y] 



Q 2 VQ 2 



(2.10) 
The interference term is the most intricate one, since various CFF combinations, which 
have different ^-dependencies, enter in the harmonics and due to the BH propagators all 
of them will contribute to the ^-integrated interference term. Utilizing the exact results, 
given in (66,67,69) and appendix B.l of [75], it can be shown that due to the (f) integration 
the transverse CFFs J r -+ disappear. Moreover, the y-dependent factor (2.10), arising from 
the BH propagators, cancel exactly in all remaining expressions and we also find a unique 
y-dependence for the net result. We also emphasize that the dominant first harmonic gets 
suppressed by 1/Q and cancels a contribution in the constant term, yielding a result that 



is proportional to x\. Finally, we add that the CFFs £ ++ and £q+ are absent in the 
unpolarized interference term. 

Let us skip here further details and quote the new results for the moderate/small-XB 
region: 



da 



BH 



^a 2 em y 2 



i-§?Hv) 



(-l)K 2 (Q 2 -t) 
t(Q 2 +t)(l-x B ) 



F Kt) ~ jh*I%(t) 



+ (D(x 



dt 



1 I r>2 



-tQ 2 (l + e 2 

da™ _ Ana 2 em y(2-y) 

dt 



Q2 V-&0--XB] 



Jut 



QUi + e 2)( 2 _ 2 y + y 2 + ^1) 1 - x B 



1 + t4-2/(1 + ^ 



KeC(J- ++ |Jb+) 



(2.11) 
(2.12) 



d<^ 



irat, 



x 



B 



[C(T ++ ,T* ++ ) + C(T- + ,T*_ + ) + e(y)C(T 0+ ,^ + )] , (2.13) 



rf* QVl + e 2 1 - ^B 

where the bi-linear C-coefficient of the (D)VCS term is given in (45) of [75] and the linear 
C-coefficient of the interference term reads 

t 



^(-^++1-^0+) 



Fi(t)n 



++ 



4M 2 



F 2 (t)£ ++ -{F 1 (t) + F 2 (t)}n 



++ 



(x B ,t,Q 2 ) + 0(x B ). 
(2.14) 

Note that here the longitudinal helicity CFFs are suppressed by an additional xb factor. 
From the equations (2.11-2.13) one immediately reads off the well-known canonical scaling 
and the characteristic y hierarchy of the BH, interference, and DVCS term, given by 



rfa E 



oc 



y 



da 1 



y_ 

dt °" Q4 



, da DVCS 1 

and oc —— 

dt Q 4 



dt ~ -tQ 2 ' 
respectively. A few further comments about the variable dependencies are in order. 

• ^-dependencies 



The power behavior in y of the BH, interference, and DVCS term is modified. More 
precisely, we have for these three terms the hierarchy 

2/(2 -y) 



y 



1+32 



' yMl-x B ) 



1 + 



2/(1 + ^ 



2 - 2y + y 2 + ^ 



2y2 J 

2~ 



1, (2.15) 



Q2 y y± i Q 2 

where both the BH and DVCS term is further separated into transverse and longitudinal 
parts. The latter is proportional to the polarization parameter s(y), which, however, in 
DVCS kinematics appears to be power suppressed. The additional ^-dependence of the 
BH term, cf. (2.10), depends on both the t/Q 2 ratio and xb- At y = and in the vicinity 
of y co \ it has the values 

1 1 11 



1 and 



i + i 



i + i 



1 _ 4^ + ^b* j 2/coi ~y 1-2/ 



(2.16) 



respectively. In the DVCS kinematics this function can be approximated by 1/(1 —y). The 

additional ^/-dependence of the interference term (2.12) is given by the rather mild concave 

function 

2-V „ 2-y 

2 - 2y + y 2 + e ^f ~ 2 - 2y + y 2 ' 
which takes the value one at both endpoints y E {0, ^1} and has a maximum of « 1.21 at 
y « 0.59. 



• ^-dependence of the BH cross section 

The kinematical factor K 2 /(—t) in the BH cross section (2.11) is proportional to (t—t m [ n )/t 
and, hence, it vanishes at the phase space boundary t —> £ m in- Thereby, the BH cross section 
(2.11) remains finite and is proportional to x^/(—t m [ n ). If we have the region — t m i n <C —t 
in mind, where — t m i n ~ x\M 2 vanishes at small xb, we will loosely say that the BH cross 
section is proportional to l/(— t). We add that in this ^-region and for xb < 0.05 the 
approximation (2.11) works on the level of one percent and better. 

• small-XB region 

At small xb the CFF behavior is governed by a possible "pomeron" exchange, which yields 
that even xb x ^(xb, £, Q 2 ) may grow for decreasing xb values. Taking the limit xb —> for 
the kinematical factors of (2.11-2.13) yields for — t m [ n <C —t the rather accurate kinematic 
expressions 



da 
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Anai 
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F 2 (t) 



AM': 



■*?(*) 



dt 

da DVCS 
dt 



-tQ 2 [y col -y}\l + (l-y) 



for y < y co i « 1 + 



t 



1, 



ira 



em 2 



Q 4 



4 {C{T ++ ,T* ++ ) + C(T- + ,F*_ + ) + e(y)C(T Q+ ,T* 0+ )] , 



(2.17) 
(2.18) 



where 



Xb C (y , J~ 



\xbV-Y 



' \x B £\ 2 + 



4M| 



xb'H 



4M| 



Ixb^I 



(x B ,t,Q 2 ) (2.19) 



with the new notation 1 



£(x B ,t,Q 2 



x B 



xb 



-S(xB,t,Q 2 



(2.20) 



x This redefinition absorbs a common prefactor xb/(2 — xs+xst/Q 2 ) of 8 that appears in all C- coefficients 
and it ensures that 8 has the same phenomenological Regge counting in the small-xs region as the other 
CFFs. It cancels the l/£ ~ (2 — xb)/xb factor that appears in the form factor in front of 8, used for the 
decomposition of the DVCS amplitude (analogously for GPD E). 



The interference term is suppressed w.r.t. DVCS cross section by an additional factor 
xb y and can be safely neglected. As one immediately realizes from these rather accurate 
kinematic approximations, the BH cross section (2.17) is kinematically enhanced at small 
—t and suppressed at small y values. However, most important is that the DVCS cross 
section (2.18) in the small-XB region grows with decreasing xb, caused by an effective 
"pomeron" exchange in the t-channel. Thus, even the relative kinematical — t(l — y)/Q 2 y 2 
suppression of the DVCS cross section w.r.t. BH one can be overcome. Moreover, the DVCS 
signal can be further experimentally enhanced by an upper y cut. However, it should be 
kept in mind that the ratio of DVCS cross section to the BH one depends on the competing 
interplay of xb, Q 2 , and t dependencies. In particular, if the DVCS cross section falls off 
much faster with increasing —t than the electromagnetic form factor i*i(£), like in the case 
of the often assumed exponential t-dependence, the ratio of DVCS cross section to BH one 
can become very small at larger — t values. 

Finally, let us quote the ^/-dependence of the ^-differential cross section (2.8) in the 
most obvious manner for general DVCS kinematics Q 2 > —t and y < y co \: 



da TOT y 2 



dt ^ t yy) dt 



y (1 - J) ^/^ + ? da™ dolf™ da-- 

dt (-> . Xi-XB)t \ ( Qi+t ,.\ i_ v+ yl + fyl dt dt yy> dt 



Q 2 +t ) y Q2+x B t 

(2.21) 
The reduced BH cross section da^ il /dt J re(y)daf u /dt for the smaller-xs region, the reduced 
interference term da™ T /dt, and the DVCS cross section da^ vcs / dt + e(y)daf j vcs / dt can be 
read off from (2.11), (2.12), and (2.13), respectively. Note that the ^-dependent factor in 
front of the interference term is given by a transverse photon flux asymmetry, 



y (1 - |) vTT^ _ c. 



-++ 



l-y+ y 4 + e ^f £__ + £++' 

Depending on the kinematics, the application of the formula (2.21) is two-fold. In the case 
that the subtraction of the BH cross section can be reliably done, the measurement of this 
subtracted cross section at three different beam energies allows in principle to separate the 
longitudinal DVCS cross section, transverse DVCS cross section, and the interference term. 
One may also utilize the y-dependence to cross-check experimentally if a BH-subtraction 
procedure is well understood. 

2.2 Relating DVCS observables to GPDs 

GPDs, denoted here generically as 

F(x, rj = €,t, fi 2 ) with F G {if, E, H, E} , 
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Figure 2. Factorization of the DVCS amplitude to leading order in perturbative QCD and to 
leading twist-two accuracy. This yields equation (2.24) that expresses CFFs (filled circle) in terms 
of GPDs (filled ellipse). 

are intricate functions that, besides depending on the partonic momentum fraction x and 
the momentum transfer squared t, depend also on the ^-channel longitudinal momentum 
fraction 77, called skewness (often denoted by £ in the literature 2 ), and on the factoriza- 
tion scale fi 2 . The unpolarized parton GPDs are called H and E [77], where the former 
(latter) GPD can be loosely associated with a proton helicity (non) conserved distribution. 
Analogous nomenclature is used for the polarized parton GPDs H and E [77]. GPDs have 
certain spectral properties [56, 78] and so their x-moments are polynomials of certain order 
in 77, with lowest moments being equal to elastic nucleon form factors. In the forward limit 
(t —¥ 0, 77 —¥ 0) H (H) reduce to the unpolarized (polarized) Parton Distribution Functions 
(PDFs), commonly called (A)q and (A)g for quarks and gluon, respectively. Furthermore, 
in the region \x\ > 77, where a parton is exchanged in the s-channel, GPDs are constrained 
by positivity conditions, which can be viewed as a consequence of a wave function overlap 
representation [79-81]. However, this GPD property is exact only to LO accuracy. To our 
best knowledge, no attempt has been undertaken to derive positivity constraints for the 
77 = x case. This implies that existing positivity constraints mostly do not apply for the 
phenomenological description of deeply virtual processes. However, as we will see below, 
they are important constraints for GPD models, e.g., as used by us in Sect. 5.2 for the 
purpose of extrapolation from the 77 = x to the 77 = case. 

DVCS observables can be exactly evaluated in terms of the helicity CFFs (2.2). To 
express them in terms of GPDs in a systematically improvable manner, it is maybe ap- 
propriate to utilize a conventionally defined GPD-inspired CFF basis, such as the one 



2 £ stands for a Bjorken-like scaling variable while 77 is a second scaling variable appearing, e.g., in doubly 
virtual Compton scattering. In deeply virtual production of photon and mesons one has 77 ~ £. Note that 
below, in Sect. 3 only, the symbol 77 will be used also to denote rapidity. 
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introduced in [72] 3 : 

T e {% £, H, £, Hs, £3, #3, 4 «t, Sr, «t, &r} • (2.22) 

Here, the CFFs H, £, %, and £ are associated with twist-two GPDs F E {i?, 22, if, ^} and 
govern the photon helicity non-flip DVCS amplitude, i.e., at leading twist-two accuracy we 
have 

J- ++ (xB,t,Q 2 )=^(x B ,t,Q 2 ) + (9(l/Q 2 ) for T^{U,£,U,£}. (2.23) 

It is ensured by the factorization theorem [63, 78, 82] that these four dominant CFFs arise 
from the convolution of twist-two GPDs with hard coefficients, which are perturbatively 
calculable as a series in the strong coupling constant a s . Presently, these coefficients are 
known to NLO accuracy in the standard minimal subtraction scheme [82-87] and to next- 
to-next-to-leading order (NNLO) accuracy in a special scheme [88, 89]. To LO they are 
calculated from the handbag diagram, depicted in Fig. 2, yielding the convolution formula 

jr {xB ^Q2^0j2j\ x __A__ T{x ^_ x} *}(*,£, t, M 2 ) for ^{~'f}, 

(2.24) 
where e^ are the fractional quark charges. The variable £ ~ #b/(2 — #b) is a conventionally 
defined Bjorken-like scaling variable, equated to the longitudinal momentum fraction in the 
t-channel, and fi 2 ~ Q 2 being the factorization scale. Note the conventional dependence 
as function of this scaling variable is in the orders of 0(1/ Q 2 ). It further reduces if one 
takes into account kinematic corrections, evaluated to twist-four accuracy at LO in a s 
[90-93]. As is well known, the ambiguity in setting the factorization scale diminishes 
in higher orders of perturbation theory as long as the perturbative corrections to the 
GPD evolution are consistently taken into account. Moreover, as long as we consider 
only the DVCS process, the perturbative order to which we describe its amplitude can be 
mainly understood as a convention (in a DVCS scheme, like in the DIS scheme, only the 
perturbatively predicted evolution would alter, if we would switch, e.g., from LO to NLO). 
In the minimal subtraction scheme the evolution kernels are known to NLO accuracy [86]. 
We also recall the fact, well known from unpolarized DIS, that the absence of gluon GPDs 
in the LO convolution equations (2.24) does not imply that these GPDs are absent from a 
LO description; they drive the evolution of the sea quarks. 

The CFFs 7/3,^3,^3, and £3 are expressed by twist-three GPDs, containing informa- 
tion on three-parton correlation functions, and enter into the photon helicity longitudinal- 



3 To simplify notation we set here F3 = 2£ (j\£* 3 — J\^ 3 ). Note that the prefactor £ does not imply 
that F3 vanishes in the £ — >• limit. 
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transversal flip amplitude, which reads to twist-three and LO in a s accuracy as 

Fo+(xB,t, Q 2 ) = y^- v [x B T + -F3] (x B ,t, Q 2 )+0 (1/Q 2 )+0(« S ) , 

QvTT^ (2 - x B + |f) 

(2.25) 

where 



x&t\ ( t N 



K = ] J-(1- xb) \l + -±)t-\l + -g) x\Ml (2.26) 

is a kinematical factor that vanishes at the minimal value of —t. The CFFs 'Ht,£t?'Hti 
and St are the dominant contributions to the transverse helicity flip DVCS amplitude, 
which, at leading twist accuracy, arises from the transversely polarized gluon GPDs that 
are perturbatively [94] and power suppressed [92, 93, 95]. Consequently, we have 

J-_ + (x B ,t,Q 2 )= ^T(x B ,t,Q 2 ) + 0(l/Q 2 ) with F T (x B ,t,Q 2 ) = 0(a s ). (2.27) 

If not stated otherwise, in the following we work for convenience to twist-two and LO 
accuracy, where we take four light quarks and we adopt the conventions 

£=-??— and i? = Q 2 . (2.28) 

2 - xb 

With these approximations GPD phenomenology can be drastically simplified. Namely, 
the convolution formula (2.24) tells us that the imaginary parts of the four dominant CFFs 
are given by the GPDs on the cross-over line x = £, 

Qm^(x Bj , t, Q 2 ) L = ttF(£, f, t, Q 2 ) , Fg {if, £7, tf , £} . (2.29) 

Furthermore, by means of the GPD spectral property one obtains from (2.24) a dispersion 
integral representation for the real parts of these CFFs [96] , 

»J H \(x Bl ,t.Q 2 ) L = PVjf'<i.T^^|^J(.T,.T,*,Q 2 )TC(t,Q 2 ), (2.30) 

Se{*~} (**, t, <?) * *vf\ ^ { «*} (x, ,, «, fi») + > { - (t ° Q2) } . (2.31) 

Here P = — C, entering in (2.30) as subtraction term, is given as convolution of the so called 
D-teim contribution (introduced in [97] to complete GPD polynomiality in one possible 
manner), which can be extracted for a given GPD. Note that the dispersion relation for S 
is over-subtracted and that the subtraction constant C(£, Q 2 ) contains a pion pole contribu- 
tion. This pole contribution can be calculated rather analogously to the D-term, e.g., from 
the suggested parameterizations [98, 99] or from extraction using a Regge-inspired GPD 
parametrization [100]. Since in this approximated framework at fixed photon virtuality 
only the GPDs at the cross-over line x = £ and two subtraction constants enter, GPD 
phenomenology is drastically simplified. 
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2.3 Present status of DVCS measurements and GPD analyzes 

Let us first consider experiments which have only an electron beam available. The three 
parts of the electroproduction cross section (2.1) contribute, depending on the kinematics, 
with different strength to the various harmonics. One can remove the BH cross section 
(2.3), taken to LO accuracy in a em , by measuring the cross section differences for single 
spin flip observables, e.g., the beam-helicity difference Alu 



dA Ll] 1 



dx B dtdQ 2 d(f) 2 



da^ da" 



(2.32) 



dx^dtdQ 2 d(f) dx^dtdQ 2 d( 
and analogously for a longitudinally (Aul) an d transversely (Aut) polarized proton target. 
These observables are expanded in terms of odd harmonics 4 . In fixed target kinematics they 
are mainly dominated by the sin(0) and/or sin(^) cos(0) harmonics of the interference term 
(2.4), giving access to the imaginary part of four twist-two associated CFF combinations, 
see (4.13-4.16) below. However, the DVCS term (2.5), suppressed in these observables by 
1/Q 2 , may also contribute to some extent. We add that in double spin flip experiments 
the BH cross section (2.3) also enters, however its cos(0) harmonic can be quite small, 
which may allow the access to the cos(0) harmonic of the interference term, i.e., three 
combinations of twist-two associated CFF combinations. 

Unpolarized electroproduction and electron-helicity dependent cross section measure- 
ments at rather large xb and small — t have been performed with small uncertainties by 
the Hall A collaboration at JLAB [17]. The measured cross section differences (2.32) is 
compatible with various GPD model predictions, see [17, 101, 102]. In the unpolarized 
case, however, the measurements at four different —t values, at Q 2 = 2.3 GeV 2 and rather 
large xb = 0.36 indicate that the DVCS cross section at these kinematics is much larger 
and drops much faster with growing —t than expected from common GPD models. As 
explained in Sect. 2.1, at small xb (large W) the "pomeron" behavior leads to the DVCS 
amplitude outgrowing the BH amplitude and as a result of the ^-integration, the interfer- 
ence term is negligibly small in this region. Therefore, at the HI [1, 3, 4, 6] and ZEUS 
[2, 5] collider experiments the DVCS cross section has been accessed by subtracting the BH 
cross section. Thereby, the subtraction method has been checked experimentally, since in 
some parts of the kinematic phase space the BH cross section dominates and Monte Carlo 
simulations can be directly confronted with measurements. The size of the cross section 
was predicted by a simple model [103] and can be at NLO also described with standard 5 



4 Here and in the following sin(n0), cos((/?) sin(n0), and also sin((/?) cos(n0) are called odd harmonics, 
while cos(n(/>), cos((^) cos(n0), and also sin((p) sin(ri(/>) are called even harmonics. 

5 We distinguish here between standard and flexible GPD models. Former, e.g., set up in [72, 104-107], 
rely on a more or less fixed skewness prescription and are used in model predictions the latter allow for a 
flexible adjustment of the skewness effect and a consistent GPD description of present DVCS data. 
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GPD models, however, not at LO [108, 109]. A simple flexible GPD model allows to de- 
scribe the HERA collider data at LO, NLO, and NNLO, which allows to quantify GPD 
reparametrization effects [110]. 

In some experiments only asymmetries, less affected by possible normalization prob- 
lems, are measurable. Having only an electron beam at hand one can access the interference 
term with single spin flip experiments by polarizing the electron beam longitudinally (elec- 
tron beam-helicity asymmetry) 

/ da^ da^ \ / / da^ da* 



\dx B dtdQ 2 d(f) dx B dtdQ 2 d(j)J / \dx B dtdQ 2 d(j) dx B dtdQ 2 dcj) J 

and analogous equations hold true for single spin flip asymmetries with longitudinally 
(Aul) or transversely (Aut) polarized nucleons and unpolarized electron beams. Here, the 
squared BH term in the numerator will drop out again at LO accuracy in a em and the 
squared DVCS term will yield some contamination, while the normalization is governed by 
all three terms of the unpolarized cross section (2.1). In addition to longitudinal proton 
spin asymmetry measurements at HERMES [10] and CLAS [14], electron beam-helicity 
asymmetries were measured at CLAS [15, 16]. 

The HERA experiments had both electrons and positrons beams available, which al- 
lowed to access the interference term via the beam charge asymmetry 

/ da + da~ \ I ( ^ a+ ^ a ~ 



A ° \dx B dtdQ 2 dc/) dx B dtdQ 2 dcf)) I \dx B dtdQ 2 dcf) + dx B dtdQ 2 dc/)) ' ^' M > 

where the numerator is entirely given by the interference term, however, the normalization 
depends also on the DVCS squared term. This asymmetry has been measured by the 
HERMES collaboration [7], where the correlation between the lowest and first harmonics, 
predicted in [72], was confirmed. The beam charge asymmetry was measured also by the 
HI collaboration [6] at large W (small x B ) where, however, this observable (as well as the 
^-differential cross section and the longitudinal spin asymmetry) is dominated by the CFF 
% and uncertainties are large. Hence, the CFF £ , giving access to sea quark and gluon 
GPD E that enters Ji's angular momentum sum rule, could not be revealed at small x. 

The HERMES collaboration provided the most complete measurement of thirty-four 
DVCS asymmetries, where a missing-mass event selection method was employed. This 
includes also a partial interference/DVCS decomposition for asymmetries measured with a 
transversely polarized [8] and unpolarized [9] proton target. However, since the normaliza- 
tion depends on the unpolarized DVCS cross section and both statistical and systematical 
uncertainties are rather large, a full disentanglement of twist-two related CFFs and an 
access to the twist-three sector could not be achieved. In particular, GPD E cannot be 
accessed from these measurements in a GPD model unbiased manner. 
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Figure 3. Results of least-squares fits in two scenarios with only a small number of CFFs locally 
fitted to data separately for each of 12 HERMES bins. First, with only ^mH and !Re£ fitted (red 
diamonds) and, second, with SmH, 3?eH and *<smH (purple pluses). For comparison, result of a 
one-to-one mapping procedure is also shown (green stars). 

This large set of DVCS observables, measured by HERMES in twelve kinematical bins 
(some are measured in 18 bins), allows for a local extraction of CFFs. Since experimental 
uncertainties are rather large for most of the observables, one may still rely on the hypoth- 
esis of twist-two dominance and extract the twist-two associated CFFs by maps [72] , by 
least-squares fits [111-114], or neural networks [115]. To avoid a misinterpretation of ex- 
perimental measurements, these local methods should be utilized with care. In particular, 
differences exist between the view points of random variable map and regression methods, 
see Fig. 3. A one-to-one map of eight twist-two dominated asymmetries into the space of 
CFFs reveals that only the imaginary part of the CFF % significantly differs from zero 
while its real part and the imaginary part of CFF % are relative small. All other twist-two 
dominated CFFs have large uncertainties and are compatible with zero [116]. By means 
of the LO approximation (2.29) the results for the imaginary parts can now be viewed as 
GPDs on the cross-over line, while the dispersion relations (2.30, 2.31) may in principle 
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be utilized as sum rules to constrain the GPDs on parts of the cross-over line that are 
outside of the accessible kinematics [117]. We add that so far no attempt has been made 
to access photon helicity flip contributions, related to twist-three and transversity GPDs. 
However, the smallness of higher harmonics is compatible with the hypothesis of twist-two 
dominance. 

Certainly, the partonic interpretation of DVCS measurements, the inclusion of the 
Q 2 evolution, perturbative corrections, kinematic corrections [90, 91], and the access to 
three-parton correlations [72, 118] requires a global analysis with flexible GPD models. 
Having measurements over a wide Q 2 range allows, through evolution, to reveal the GPD 
away from the cross-over line. This is used for the description of the DVCS cross section 
measurements at small xb, whereas for fixed target kinematics the Q 2 lever arm is small 
and evolution effects are relatively weak (for an example study see [119]). 

In a first step of a global DVCS analysis, unpolarized proton data were employed in 
GPD fits [110, 119, 120]. Thereby, the world data set could be described with x 2 /d.o.f. « 
1, using the KM10 model. Nevertheless, in such a fit the four CFFs 1~L,£,T-L,£ cannot 
be disentangled and, partially for this reason, even the dominant T~L suffers from larger 
uncertainties, see Fig. 4. Below we will also employ the model KMlOa, which has also a 
good x 2 /d.o.f. « 1 fit to the data set, but ignores the Hall A cross section measurements. 
Including polarized proton data in a global fit could certainly help to disentangle CFFs 
even better. In a more recent fit, given in [116], we found that the KM model, designed 
for the unpolarized case, describes even such a set of DVCS data with x 2 /d.o.f. « 1.6, 
where most of the tension is due to the four unpolarized cross section measurements of 
Hall A collaboration. We emphasize that this tension can have different origins, e.g., it is 
maybe prudent to still consider the possibility that the experimental issue of exclusivity 
plays a role in most of the world data. For instance, beam spin asymmetry measurements 
from the HERMES collaboration with a complete event reconstruction yields an increase 
of their size, softening, thereby, the tension between measurements and standard GPD 
predictions [12, 102]. ^From present DVCS data, we can certainly state that GPD H plays 
the dominant role, some phenomenological constraints for GPD H can be obtained, and 
proton helicity flip GPDs E and E remain unconstrained. Let us add that present GPD 
phenomenology includes also deeply virtual meson production, in first place in the hand- 
bag model approach [106, 107, 121, 122] and was started in the perturbative factorization 
framework with flexible GPD models [100, 123]. So far a reasonable description of the 
considered deeply virtual meson production channels and DVCS, currently explored on the 
level of LO accuracy, can be reached [102, 119, 123] except for the large-^B region. 

As pointed out and illustrated in Figs. 3 and 4, present DVCS measurements provide 
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Figure 4. ^smH/ir obtained from DVCS observables with different strategies: hybrid model fits 
KM 10 (solid) KMlOa (dashed), KMlOb (dash-dotted) [Hall A cross section data are neglected], 
GK07 model from DVEM (dotted) [106], seven-fold CFF fit [111, 112] with boundary conditions 
(squares), H : % CFF fit [113] (diamonds), smeared conformal partial wave model fit [124] within 
H GPD (circles). The triangles result from neural network fit [115]. 

some limited information on GPDs and future precision measurements are required to pin 
them down. New fixed target experiments are planned at COMPASS-II with a polarized 
muon beam, extending the HERMES kinematics to lower xb, and JLAB-12 GeV will 
bridge the gap between the kinematics of the present JLAB experiments' to the HERMES 
experiment, see Fig. 5. Moreover, a high luminosity machine in the collider mode with 
polarized electron and proton or ion beams has been proposed [73] and will be introduced 
in the next Section. 

3 The EIC project and Monte Carlo simulation 

In order to open a new window into a kinematic regime that allows the systematic study 
of quarks and gluons, EIC is designed to provide a wide range in c.o.m. energies, polarized 
lepton and light ions beams and heavy ion beams, all at a very high luminosity [73]. This 
creates an unprecedented opportunity for discovery and precision measurements, and would 
allow us to study the momentum and space-time distribution of gluons and sea quarks in 
nucleons and nuclei [73, 125]. The main requirements for an EIC machine are: 

• Highly polarized (> 70%) electron and proton/light ion beams; 

• Ion beams from deuteron to heaviest nuclei (uranium, lead); 

• Variable center of mass energy, ranging from about 20 GeV up to 150 GeV; 

• Collision luminosity ~ 10 33-34 cm _2 s _1 . 
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Figure 5. The kinematic reach in x — Q 2 for existing DVCS measurements from HI /ZEUS, 
HERMES, CLAS and Hall A, as well as planned ones at COMPASS II and JLAB@12GeV, and the 
proposed EIC. 





Figure 6. The Layout of the two proposed EIC machines: eRHIC (left) and ELIC (right). 

Two independent designs for a future EIC have evolved, eRHIC and ELIC, both using 
part of already available infrastructure and facilities (see chapter 5 in [73]). At Brookhaven 
National Laboratory (BNL) the eRHIC design (Figure 6 left) utilizes a new electron beam 
facility based on an Energy Recovery LINAC (ERL) to be built inside the RHIC tunnel to 
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collide with RHICs high-energy polarized proton and nuclear beams. At JLAB the ELIC 
design (Figure 6 right) employs a new electron and ion collider ring complex together 
with the 12GeV upgraded CEBAF, now under construction, to achieve similar collision 
parameters. The kinematic phase space achievable at an EIC for electron-proton collisions 
is shown in Fig. 5 and compared to existing DVCS data and planned future experiments. 
At an EIC it will be possible to study DVCS measuring, for the first time simultaneously 
and with high accuracy, both differential cross section and spin and charge asymmetries in 
a kinematic range that extends from large xb, typical for fixed target experiments, down 
to small xb, typical for the HERA collider experiments. 

The present study is based on the eRHIC version of an EIC and its new dedicated 
detector, designed to fulfill the requirements for the golden experiments at an EIC and thus 
being simultaneously highly efficient for inclusive, semi-inclusive and exclusive reactions. 
The eRHIC expected luminosity for ep collisions as a function of the beam-energy is shown 
in Figure 7. At eRHIC the full range of proton-beam energies will be at hand from the 
early beginning of operations, whereas the energy of the new electron-beam will be initially 
at 5 — 10 GeV (stage I) and will be later upgraded to higher energies up to 20 — 30 GeV 
(stage II). The newly designed eRHIC detector, shown in Fig. 7, will have the following 
properties: 



3T Solenoid 




E p (GeV) 



hadron-beam 



Figure 7. Left: the expected luminosity at the eRHIC collider as a function of the beam-energy 
configuration for ep collisions. Right: a sketch of the eRHIC detector. 

• Wide acceptance — 5 < 77 < 5 for both the scattered lepton and the produced hadrons; 

• The same rapidity coverage in electromagnetic calorimetry and tracking; 

• High electron track finding/reconstruction efficiency, capability to discriminate two 
electromagnetic clusters down to a difference of 1 degree of polar angle in the rear 
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endcap electromagnetic calorimeter and good precision for momentum (energy) re- 
construction; 

• Particle identification to separate electrons and hadrons as well as pions, kaons and 
protons over a momentum range of 0.5 GeV to lOGeV for rapidities between -1 to 1 
and 0.5 GeV to 80 GeV for 1< \r)\ < 3; 

• Good vertex resolution; 

• High acceptance for forward going protons and neutrons from exclusive reactions as 
well as from heavy ion breakup (Roman Pots and Zero Degree Calorimeter will be 
part of the detector). 

• Low material budget to reduce electron bremsstrahlung and to achieve good resolu- 
tion in the reconstruction of all the kinematic variables. 

• Very small low scattering angle forward scattered electron tagger (Q 2 < 0.1 GeV 2 ) 

The Monte Carlo (MC) generator used in the present study is MILOU [126], which 
simulates both the DVCS and the BH processes together with their interference term. It 
is based on the code by Freund/McDermott [127, 128], which utilizes the approximations 
described in [72], and is tuned to HI and ZEUS measurements. The DVCS amplitude is 
evaluated in a GPD-inspired framework to NLO accuracy [82-85] , including the NLO GPD 
evolution [86], by a routine, which provides tables of CFFs. The real and imaginary parts 
of CFFs then are used to calculate the cross sections for DVCS, BH and their interference 
term. The t-dependence of the DVCS amplitude is introduced as an exponential, i.e., the 
DVCS cross section reads 

da DVCS (W^,Q 2 ) 



oc 



exp{B(Q 2 )t} 



dt 

with the exponential t-slope parameter B(Q 2 ) being constant or having a logarithmic Q 2 - 
dependence 6 . The MILOU code has been slightly modified from its original version as 
described in Appendix A. 

The simulations used for our studies are based on the following MILOU options: 

• The slope B(Q 2 ) = 5.6 GeV -2 is set to be constant. 

• CFFs tables are generated from a GPD model to NLO and twist-two accuracy. 



6 If the t-dependence of flavor singlet quark and gluon GPDs is chosen differently at the input scale, 
perturbative evolution will alter the t-dependence for the resulting DVCS cross section. This should not 
be confused with the MILOU option to alter additionally the Q 2 -dependence of the exponential t-slope by 
hand for a given GPD model. 
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• Proton dissociation background, ep — >► e^jY ', has not been included in the simulation. 

To our best knowledge, the first two choices guarantee that a pure and consistent GPD 
framework is utilized in the MILOU simulations, see [129] and footnote 6. 

The DVCS and BH processes have been simulated according to the following selection 
criteria: 

• Q 2 > 1 GeV 2 ; 1CT 5 < xb < 10 _1 ; binned logarithmically in 4 Q 2 - and 5 xe-bins per 
decade and in several |i|-bins; the bins in Q 2 are: 1.0 < Q 2 < 1.78 GeV 2 ; 1.78 < Q 2 < 
3.16 GeV 2 ; 3.16 < Q 2 < 5.62 GeV 2 ; 5.62 < Q 2 < 10 GeV 2 ; 10 < Q 2 < 17.78 GeV 2 . 

• Detector acceptance criteria: 0.01 < y < 0.85 and \rj\ < 5.0 for scattered electron 
and produced photon, scattered proton acceptance: 0.03 < |t| < 1.5 GeV 2 (proton 
detected in the roman pots); 

• BH rejection criteria applied for the cross section measurement: em-clusters-energy 
> lGeV; 9 d -0 1 >0rad. 

The Q 2 and xb range is within the phase space reachable with an EIC/eRHIC. The electron 
and proton beam-energy configuration considered for the present study are: 5 x 100 GeV 2 , 
5 x 250 GeV 2 (for stage I) and 20 x 250 GeV 2 (an example for stage II). 
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Figure 8. The distribution of statistics in each {Q 2 ,£b} bin for eRHIC stage I (left and middle) 
and stage II (right) at a luminosity of 10fb _1 . 

For the purpose of DVCS cross section measurements it is important to remove from 
the signal the background coming from the BH events. The latter is a QED process, well 
known to an uncertainty of the order of 3% coming from the uncertainty on the proton 
form factors. It can be subtracted from the signal by means of a MC technique. Thus, 
especially at a high luminosity machine like eRHIC where systematic uncertainties will 
dominate the measurements, it is important to minimize the BH contribution, particularly 
at low c.o.m. energy, where BH tends to dominate over the DVCS (see Sect. 2.1). The 
fraction of BH events has been estimated using a MC sample containing both DVCS and 
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BH processes. The BH contamination was investigated for each {Q 2 , xb, i] bin as a function 
of the electron energy loss y. After all BH suppression criteria have been applied it was 
found that at large c.o.m. energies the BH contamination grows from negligible (at low-?/) to 
about 70% at y ~ 0.6 allowing for a safe BH subtraction, whereas for lower c.o.m. energies 
the BH contamination grows faster with y and can be dominant depending on the bin; 
nevertheless most of the statistics at this low c.o.m. energy is contained in the safe region 
y < 0.3. 

Figure 8 compares the distribution of the statistics per bin for the eRHIC beam-energy 
configurations 5 x 100 GeV 2 , 5 x 250 GeV 2 (both reachable at a stage I) and 20 x 250 GeV 2 
(available at a stage II), considering an integrated luminosity of ~ 10fb _1 . The results 
shown in the present paper are based on simulated data samples corresponding to an 
integrated luminosity of 100fb _1 for the 20 x 250 GeV 2 configuration and 10fb _1 for the 
5 x 100 GeV 2 configuration, both corresponding to approximately 1 year of data taking at 
eRHIC assuming a 50% operational efficiency. The data samples generated for the propose 
of measuring the differential cross section only contain the DVCS process whereas samples 
containing DVCS, BH, and their interference term have been generated for measurements 
of different single spin asymmetries. 

All the generated events have been smeared according to expected momentum and 
angular resolutions. The statistical uncertainty for the differential cross section can be 
at small values of —t as low as few percent; the same is true for the uncertainty for the 
extracted slope parameter B. This implies that the measurement is actually limited by 
systematics. For the purposes of the present work, a systematic uncertainty of 5% has been 
assumed, based on the experience at HERA and the expected coverage and technology 
improvements of the new detector at eRHIC. The overall systematic uncertainty, due to 
the uncertainty on the measurement of luminosity, is not considered for this paper as it 
simply affects the normalization of the cross section measurement. 

4 Selected DVCS observables at EIC 

As explained in Sect. 2, the isolation of CFFs is a rather intricate task, which can be 
only achieved by measuring a complete set of observables. However, we have also seen 
that photon helicity flip contributions, which are suppressed in DVCS kinematics, are not 
traceable in the present world data set. Hence, we restrict ourselves to four twist-two 
associated CFFs to study the physics case of DVCS measurements at a suggested eRHIC, 
giving emphasis to twist-two dominated observables. As motivated in Sect. 3, we choose two 
scenarios: one with a relatively low and another with a high c.o.m. energy, corresponding 
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to the beam configurations 

E e xE p = 5x lOOGeV 2 and E e x E p = 20 x 250GeV 2 . 

For future DVCS measurements at 5 x 100 GeV 2 it is maybe expected that the de- 
scription of precise data in this region of transition to the small-xs physics requires rather 
complex GPD models, which are not needed for the description of the present DVCS data. 
For the higher energy case it is expected that valence quark contributions are negligibly 
small and non-negligible CFFs 

xb x % and potentially xb x £ 

are governed by an effective "pomeron" exchange in the ^-channel, associated with both 
sea quarks and gluon contributions, and that they (moderately) grow with decreasing xb- 
Thereby, almost nothing is known about the CFF £ , which, as pointed out in Sect. 2.3, is 
not accessible from present DVCS measurements in neither the collider nor the fixed target 
mode. The available theoretical/phenomenological guidance is not yet fully trustworthy. 
On one hand a "pomeron" coupling to proton helicity non-conserved quantities such as the 
CFF £ is phenomenologically not established, see Ref. [130] and references therein. On 
the other hand a pomeron like behavior for the CFF £ is perturbatively predicted by GPD 
evolution 7 . A separate study on the access of GPD E as well as the transverse spatial 
distribution of sea quarks and gluons at stage II will be presented in Sect. 5, which without 
additional information or assumptions is hard to achieve for EIC measurements at low 
beam energies. 

We expect that the remaining two twist-two associated CFFs if multiplied with xb, 

xb x T-L and xb x £ « xb x £ , 

2 - x B 

go to zero in the small-^B region. Note, however, that in contrast to the CFF H, Regge 
phenomenology provides no clear guidance for their small-XB behavior. The phenomenolog- 
ical situation is analogous to the polarized DIS function g\ (GPD H embeds the polarized 
PDF Ag). We emphasize that these essentially unknown contributions may play a role at 
the stage I kinematics. 

To cover possible scenarios, we employ in our studies three hybrid models, where sea 
quark and gluonic components of CFFs 1~L and £ are based on GPD models that include 



7 As for the perturbative evolution of unpolarized PDFs in the flavor singlet sector, the evolution of both 
GPD H and E in this sector is at small x driven by gluons, which generate an effective 'pomeron' like 
behavior. The solution of the evolution equation yields in fact an essential singularity rather a pole. Such 
a behavior can be only avoided if both the quark singlet and gluon GPDs vanish simultaneously. 
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the perturbative evolution, while their valence quarks and remaining GPDs are treated 
with dispersion relations as described in Sect. 2.2 . Two of the models are pinned down 
from global fits to the world data of unpolarized DVCS measurements, which are described 
very well, despite having rather different partonic content. We now list the models and 
describe their main properties. 

• KM 10 describes the world data set of DVCS measurements using an unpolarized 
proton target. It contains the twist-two GPDs H and H, while the real part of 
helicity-flip CFFs £ and £ are only given by subtraction constants in the dispersion 
relation (related to so-called D-term and pion pole contribution, respectively). Both 
GPD H and the (real) CFF £ are rather large and they are considered as effective 
degrees of freedom that allow to describe the unpolarized cross section measurements 
from the Hall A collaboration [17]. 

• KM 10a is analogous to the KM 10 model; however, the Hall A cross section measure- 
ments are not well described. In this model the GPD H is the dominant one, H is 
set to zero, and £ contains only the pion pole, which is accounted in the standard 
way [98, 99]. 

• AFKM12 is a flexible GPD model for the small- x region, specifically designed for the 
present study. It contains besides the sea quark and gluon GPDs i7 sea and H G also 
a flexible small- x parametrization of GPDs £" sea and E G . All of these GPDs include 
a "pomeron" behavior, which can be individually adjusted at the input scale. The 
normalization of E-type GPDs is controlled by the anomalous magnetic moment of 
sea quarks ft sea = 1.5, which is fixed to be positive and rather large. The parton 
polarized GPDs H and E are set to zero. 

Our small- £ GPD models are set up in terms of (conformal) GPD moments rather 
than in x-space, at the input scale Q 2 = 4GeV 2 for four light quarks. They yield, similarly 
to other GPD models, the following effective functional form 8 of CFFs: 

'na(t, Q 2 



n , , . o2 



£ f (xB,t,Q?)~ir 



cot 



r a(t ' Q2) \ ha } (*, (?) , (4.i) 



2 
which resembles a Regge phenomenological ansatz with a linear "pomeron" trajectory 

a(t) = a(t = 0) + Oct. (4.2) 



8 This form arises exactly in the small- x limit of standard GPD models at the input scale; however, 
strictly spoken it is not stable under perturbative evolution. Nevertheless, the resulting CFF output of a 
GPD model can be reparameterized for a given Q 2 value and put in a Regge- inspired form. Thereby, the 
"pomeron" trajectory is altered, indicated by its Q 2 -dependence. 
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Table 1. Some selected model parameters for unpolarized sea quark and gluon GPDs (upper 
table), valence H and H GPDs as well as for subtraction constants (lower table), where squared 
mass parameters are given in GeV 2 and slope parameters a' and B in GeV -2 . 



In the KM10 and KM 10a models a dipole parametrization (l — -^) for the residual 
^-dependency was taken, while the AFKM12 model alternatively relies, as in the MILOU 
simulation, on an exponential ansatz e bt . The boundary value of the residue h a at t = 
depends on both the momentum fractions N l (Q 2 ), carried by the unpolarized parton type 
z, and the skewness effect, parameterized in terms of two model parameters s\ and s\, which 
control both the normalization of the CFFs and their Q 2 evolution; a detailed discussion 
is given in [110, 131, 132]. Analogously, we parameterize in the AFKM12 model the GPD 
E with an independent set of parameters, however, here the momentum fractions A^(Q 2 ) 
are replaced by the partonic gravitomagnetic moments B l — N 1 k % , parameterized at the 
input scale by the product of the momentum fraction N l and the anomalous magnetic 
moments n l . The momentum and gravitomagnetic sum rules are utilized to fix the gluonic 
momentum fraction and gravitomagnetic moment, respectively. From a DIS fit the PDF- 
related parameters were found [110], 



N s 



0.152. 



a 



af a = 1.158, 



a 



G 



a 



G 



g = 1.247, 



(4.3) 



which we, for simplicity, also adopt for GPD E in the AFKM12 model. Some other relevant 
model parameters are listed in Tab. 1, where, again for simplicity, we equate the Regge 
slope parameters of GPD H and the residue slope parameter for GPD E with those of 
GPD if, 



a 



i G 



a 



I sea 



Lsea 
°E 



°E 



(4.4) 



Finally, we specify the remaining GPDs on the cross-over line and the form of sub- 
traction constants, where the CFFs are calculated from (2.29, 2.30, 2.31). Only the target 
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helicity conserved GPDs on the cross-over line are modeled 

rva i/ .x 1.35 r / 2x \ _Q: ^/1 — x\ /„ 1 — x 



fl-.(», x, = — ^— j ^— j ^1 - — ^ j , (4.5) 

, 0.6F / 2x \" Q(i) (l-x\ b ( 1-xtY 1 tA s 

Here, the skewness effect is parameterized by the ratios 

v H(x,x,0) _ _ .. H(x,x,0) 

r = lim — — and r = lim — , 

x^o if (a, 0,0) ^Off(x,0,0) 

a(t) = 0.43 + 0.85£/GeV 2 , b (b) controls the x -> 1 limit, M val ( M) the residual t- 
dependence, where q(x) = H(x, 0, 0) (Ag(x) = i7(x,0,0)) are unpolarized (polarized) 
reference PDFs, e.g., the LO parametrization of [134] ([135]). The subtraction constant is 
normalized by c (rv) and the cut-off mass M c (M^) controls the t-dependence: 

C(t) - , C , 2 , C(t) = 2A6 * r * x2 , (4.7) 

!-lfe) ("4-*)( 1 "]0? 

where m^ « 0.14 GeV is the pion mass and the normalization factor 2.164 in the pion pole 
contribution matches the residue of the t = m 2 pole from the pseudo scalar form factor 
2g A M*/(l + ml/ Ml) with M n = 1.17 GeV. Note, however, that in the GPD framework 
the normalization of the pion pole contribution remains unknown. In the KMlOa model 
we use the pion pole parametrization of [133]. More explanations on these simple param- 
eterizations can be found in [110]. The parameters of the KM10 and KMlOa models are 
listed in Tab. 1. 

In the remainder we illuminate the richness of a possible experimental DVCS program 
at an suggested EIC. Thereby, we will concentrate on observables that are dominated by 
twist-two associated CFFs. In Sect. 4.1 we restrict ourselves to the unpolarized cross section 
and in Sect. 4.2 to single spin asymmetry measurements. In Sect. 4.3 we will comment on 
further DVCS related measurements, which are interesting on their own, and we shortly 
discuss the use of an unpolarized positron beam to disentangle photon helicity non-flip 
contributions from longitudinal- transverse helicity ones. 

4.1 Cross section measurements at stage I 

As emphasized in Sect. 2.1, the separation of the measurable electroproduction cross section 
(2.8) into its three parts in the most model independent way and/or with a minimal set 
of assumptions is an important goal. So far the extraction of the ^-differential DVCS cross 
section, entering in (2.8), has been only reached in the small-XB and 0.1 GeV 2 < — t < 
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0.8 GeV 2 region by the HI and ZEUS collaborations. Thereby, the subtraction method 

rfq DVCS (x B ^Q 2 ) ^ da^(x B ^Q 2 ) _ da^(x B ^Q 2 ) 

dt ~ dt dt ^ ' ' 

was utilized, where the interference term could be safely neglected and the BH cross section 
was simulated. The latter was cross-checked experimentally in the BH dominated phase 
space region. 

To understand whether such a subtraction procedure would be also reliable in the 
EIC kinematics and whether one can improve this method by utilizing the variable beam 
energy option, we consider first the generic dependence of the ^-differential cross section 
(2.8) on its variables. According to what was explained in Sect. 2.1, for smaller value of 
—t S> — tmin ~ x b^p an< ^ l ar g e V the BH cross section dominates, since it is enhanced 
by the kinematical prefactor y 2 /(—tQ 2 ). On the other hand in the limit y — >► both the 
BH cross section and the interference term drop out, where e(y — 0) = 1 and, thus, the 
sum of the transverse and longitudinal DVCS cross sections can be accessed, see (2.21). 
Moreover, the interference term (2.12) has the same canonical 1/Q 4 scaling as the DVCS 
cross section (2.13), however, it has an additional prefactor x B y. Restricting ourselves to 
the dominant CFF H, we find that the ratio of interference term (2.12, 2.14) to the sum 
of BH (2.17) and DVCS (2.18) cross sections is estimated, for smaller-xe values, to be 



da™ / _(i _ y)t F L (t) ^z# Mex B n(x B , t, Q 2 ) 

da ™ + da ovcs - 2 ^V q2 m _ _^ m + __^ {xBn{x ^ t Q2)|2 • (4-9) 

Obviously, the suppression factor 2xb \/— (1 — y)t/Q 2 < xb (DVCS requires — t <C Q 2 ) 
makes this ratio small for EIC kinematics. Moreover, we expect from Regge arguments, 
consistent with phenomenological findings, that the real part of the dominant CFF % is in 
the small- and even moderate-XB region much smaller than its imaginary part (at least for 
smaller values of — £, see the results from HERMES in Fig. 3). We conclude that in most 
of the stage I bins, given in Sect. 3, the interference term is negligible and we can simplify 
the ^-differential cross section (2.21) to 



v< 2 
da 



TOT y 



7 BH,red , BH,red 

dt ^ t \y) dt 



+ ^ (4-10) 



dt (i-W) (££*-*) " dt 

da DVCS( \ da DVCS d(J DVCS 

with dt yy) — — ^ \- e(y) — ^ — . The smallness of the interference term has been 

also seen in numerical GPD model calculations. Thereby, the use of the approximate equa- 
tions in [72] naturally yields only incomplete cancelations in the ^-integrated interference 
term. This causes the ratio (4.9) to appear proportional to {—t/Q 2 )^' 2 rather than to 
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xb x (— t/Q 2 ) 1 / 2 . Nevertheless, also in the MILOU simulations, based on the approximate 
equations in [72], the interference term turns out to be negligibly small. 

For an EIC experiment the equation (4.10) provides a further handle to cross-check 
experimentally the BH subtraction procedure. However, we expect that a Rosenbluth 
separation of the transverse and longitudinal DVCS cross section will be difficult to achieve 
in the small —t region. To suppress the BH contribution a relatively small y is needed, 
which also means that the variation of e(y), which functional dependence can be mimicked 

2 3 

by a truncated Taylor expansion e(y) ~ 1 — ^ — ^-, is only small. Moreover, if we stick to 
the twist-two expansion of the DVCS amplitude, the longitudinal DVCS cross section in 
the small-^B region will be expressed by twist-three associated CFFs and this cross section 
will be kinematical suppressed by a factor K 2 /Q 2 « — t/Q 2 , see (2.13), (2.25) and (2.25). 
On the other hand these behaviors may offer the possibility of access to the twist-three 
contribution at larger values of — t, which, in turn, allows the variation of y over a larger 
region. However, such an access may only be possible if the ^-dependence of CFFs, as 
compared to that of electromagnetic form factors, is rather flat. 

The transverse DVCS cross section contains both non-flip and transverse flip helicity 
amplitudes, where the latter would be perturbatively suppressed by (a s /27v) 2 or 1/Q 2 
corrections. Neglecting the suppressed photon helicity flip contributions and switching to 
the GPD-inspired CFF basis (2.22), we can approximately write the DVCS cross section 
for stage I kinematics as 



da BVCS _ 7ra 2 em 
dt ~ Q 4 



\xbH\ t\x-&£\ \x-qH\ t\x&£\ XB^Rel-lS* 



(1_EB) 2 AM 2 (1-^b) 2 4M 2 1-3* 



L p 2 



{xB,t,Q 2 ) 



where the functional form arises from the exact C-coefflcients by neglecting kinematically 
suppressed contributions of order O(x^) and 0{x^t/Q 2 ). As somehow expected, in our 
numerical studies it turned out that the DVCS cross section (4.11) for the 5 x 100 GeV 2 
beam configuration is rather sensitive to the choice of GPD model. To some extent this 
is also true for higher c.o.m. energies in the large —t region. In other words, a definite 
conclusion whether the subtraction method in these kinematics will be possible, cannot be 
taken without actual data. 

As mentioned in Sect. 3, the eRHIC option allows also at stage I to increase the proton 
beam energy, for kinematical coverage see Fig. 8. To illustrate the energy dependence of 
the DVCS cross section, we consider its ratio to the measurable electroproduction cross 
section 

d(7 Dvcs _ j_^d4> d(T dxBd tdQ*d4> (4 n) 
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Figure 9. Cross section ratio (4.11) of the DVCS cross section to the photon electroproduction cross 
section (2.21) as function of —t for two x B — Q 2 -bins and two different beam energy combinations. 

Considering again the CFF % as the dominant one and sticking to the small-XB approxi- 
mation with —t ^> t m [ n and y < y co \ « 1, we can estimate this ratio as 



da z 



-*(i-y) 



'i^\xBn(x B ^Q' 



da TOT 



F ?(t) ~ ^m) + =$cfi I*b«(*b, t, Q 2 )| 2 



(4.12) 



Clearly, as long as we stay away from — £ m i n , which is at EIC not reachable in the considered 
bins, this ratio will get very small at low —t and its behavior at large —t < Q 2 depends 
on the — t drop-off of CFFs. In Fig. 9 we show the typical t-shape of this ratio for an 
exponential ^-dependence at 5 x 100 GeV 2 (circles) and 5 x 250 GeV 2 (triangles) for two 
{^B, Q 2 }-bins, 



6.31 
3.98 



10" 3 <x B < 1.00 
10" 3 < x B < 6.31 



10" 2 and 3.16 GeV 2 < Q 2 < 5.62 GeV 2 (left) , 
10" 3 and 5.62 GeV 2 < Q 2 < 10.0 GeV 2 (right) . 



These results were simulated by MILOU, as described in Sect. 3. The statistical uncertain- 
ties are obtained including all the selection criteria to suppress the BH cross section also in 
the region where the DVCS cross section is extremely small. Clearly, the functional multi- 
variable dependencies, that are expected from the approximation (4.12), can be easily seen 
in the plots. In this specific GPD model, utilized in MILOU, the DVCS cross section is 
only accessible in a smaller set of {x B , Q 2 , £}-bins. However, as is clearly illustrated in Fig. 
9, an increase of the proton beam energy from 100 GeV to 250 GeV allows to overcome 
such a potential limitation. 

We take now the KM10 and AFKM12 predictions to illustrate that the DVCS cross 
section can be possibly obtained by a subtraction procedure (4.8) even at the low beam 
energy configuration 5 x 100 GeV 2 . Generally, these DVCS cross section predictions over- 
shoot those of the MILOU simulations, on the other hand the KMlOa model predictions 
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Figure 10. KM10 (left panel) and AFKM12 (right panel) model predictions of the differential 
DVCS cross section versus —t for unpolarized beams with energies E e x E p = 5 x 100 GeV and a 
luminosity of 10 fb~ . The uncertainties of the EIC pseudo data contain statistical, 5% systematical, 
and uncertainties due to BH cross section subtraction, where for the latter a 3% uncertainty of the 
BH cross section has been assumed. 

are in agreement 9 . Based on the MILOU simulations, described in Sect. 3, we obtain the 
statistical uncertainties for the model predictions by rescaling according to the ratios of 
the DVCS cross sections. All uncertainties (statistical, 5% systematical, and subtraction 
uncertainty from a 3% error of the BH cross section) were added in quadrature and the 
predicted cross section for a kinematical point, given by the center of a three dimensional 
{%B, Q 2 ,t}-bin, was assumed to be normally distributed. 

In Fig. 10 we show the KM10 (left panel) and AFKM12 (right panel) model predictions 
for the DVCS cross section versus —t for the 5 x 100 GeV 2 beam energy configuration 
for four Q 2 and three xb bins. Apart from the different t-behavior, one notices model 
differences in the normalization at lower —t values, in particular for the largest xb values. 
One also realizes that in the KM 10 model the cross section does not necessarily grow with 
decreasing xb as it is the case for AFKM12 model (solid curves), containing only the sea 
quarks and gluon components of the CFFs % and £ . Both of these observations indicate 
that valence-like contributions to 1~L and/or non-dominant CFFs can play a certain role at 



9 In the majority of bins the KM 10a and MILOU cross sections are comparable to each other, while in 
some low Q 2 and large-xs bins the KM 10a model prediction overshoots the MILOU prediction up to 100%, 
which could be attributed to model differences. 
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lower c.o.m. energies. In both panels the sizable uncertainties arise from the uncertainty of 
the BH cross section and, as expected, they appear for the low xb bins, essentially, in the 
small —t region and large —t region. Note in Fig. 10 bins are not shown in which the DVCS 
cross section is entirely dominated by the subtraction uncertainties, i.e., we ignored bins 
with y > 0.25- • -0.4. For values —t > 0.8 GeV 2 (not shown) the uncertainties associated 
with the BH subtraction become also large, particularly for AFKM12 model which possesses 
an exponential ^-dependence. We remind that all models, including MILOU, describe the 
Hl/ZEUS DVCS cross sections measurements very well (see left panel of Fig. 13) for which 
the aforementioned contributions play a minor role. 

Let us summarize the lessons for an unpolarized DVCS cross section measurement 
at rather low EIC energies. Certainly, it is safe to expect that the electroproduction 
cross sections, i.e., containing all three terms, are large enough to provide precise data, at 
present not available in this kinematical region of transition to small xb- Such data can 
be immediately included in global GPD fits; however, model assumptions will affect the 
partonic interpretation of such measurements. The isolation of the DVCS cross section is 
probably only feasible in a limited phase space (lower y values, limited —t values). Even 
in the case that this problem can be overcome by a (partial) Rosenbluth separation, the 
measurements would only provide a very qualitative insight in the transverse distribution 
of partons, since the separation of different CFF contributions is based on assumptions. 
Hence, a measurement of further observables is needed, which allows for a separation of 
the various CFFs contributions. 

4.2 Single spin asymmetry measurements 

Measuring the differences of spin-dependent cross sections (2.32) for unpolarized, longi- 
tudinally and transversely polarized protons allows the access of the imaginary part of 
CFFs in a much cleaner manner than utilizing asymmetries. In such measurements one 
may use harmonic analysis to access the imaginary parts of twist-two associated CFFs 
from the first odd harmonics, occurring from the interference of the BH and DVCS ampli- 
tudes. However, even these observables are contaminated by power-suppressed helicity flip 
contributions that stem from both the interference and DVCS squared term. The latter 
contamination can be eliminated if lepton beams of both charges are available, see discus- 
sion in the next section. This allows then for a harmonic analysis, aiming to isolate the 
imaginary parts of twist-two associated CFFs from the remaining ones. In this way one can 
separate to some extent twist-two, twist-three, and gluon transversity contributions. What 
is the best strategy to analyze a high quality data set, measured in an experiment where 
only an electron beam is available, is not so obvious at present. One may hope that, as in 
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the case of unpolarized electroproduction cross section, considered in Sect. 2.1, a common 
Fourier analysis will finally yield some simplifications and may even allow to employ the 
Rosenbluth separation method to some extent. 

For purpose of illustration we focus in the following on twist-two GPD model predic- 
tions for single spin asymmetries rather than on spin-dependent cross section differences 
(2.32). In Fig. 11 we show pseudo data that are generated using the KM10 model, and 
randomized according to the uncertainties as specified in Sect. 3 (rescaled statistical errors 
from MILOU simulations, 5% systematical uncertainty on cross section level, 5% normal- 
ization uncertainty for the polarization measurement). The error propagation from the 
^-dependent cross section to harmonic amplitudes was simply done by fitting. Note that 
the uncertainty for the projection asymptotically scales for the N 0-bins as 1/vTV, ex- 
cept for the zeroth harmonic for which scaling is l/\/2N . We note that the polarization 
error should be treated as an overall normalization uncertainty, which, however, was not 
done here. Hence, the projections on the first harmonic in Fig. 11 have an additional 
normalization uncertainty, essentially given by the polarization uncertainty. 

The upper panels in Fig. 11 show for a proton beam the electron beam spin asymmetry 
(2.33) as function of the azimuthal angle (j) for one selected bin with 5 x 100 GeV 2 beam 
energies (left panel), its projection on the dominant first sin0 harmonic, 



2-2y + y 2 Vy 2 Q 2 



^sm0 y v - 1 - y / b „ c^- 
A hv oc ~ — ~ — ; — ^ \ I i^9 x X B ^m 



FiH-^F 2 £+°f(F 1+ F 2 )H 



(x B ,t,Q 2 ) + --- , 

(4.13) 
as function of —t (middle panel), and versus xb for a low Q 2 = 2.5 GeV 2 and a high 
Q 2 = 13.9 GeV 2 value (right panel). The asymmetry is dominated by helicity conserved 
CFF % and proportional to the electron energy loss y. Consequently, if y is not too 
small, the asymmetry might be rather sizable over a large kinematical region, shown for 
5 x 100 GeV 2 (squares, thick curves) and 20 x 250 GeV 2 (triangles, thin curves). The CFF £ 
appears with a kinematic suppression factor t/4M 2 , induced by a proton helicity flip, and 
remaining CFFs also contribute, which is in (4.13) indicated by the ellipsis that include 
also further kinematically suppressed contributions. Comparing the different predictions of 
the KM10 (dashed), KMlOa (dotted), and AFKM12 (solid) models, one realizes that the 
contaminations of this asymmetry by other CFFs are in fact small. Our £ enhanced model 
prediction only slightly differs from the other ones at t ~ 0.5 GeV 2 . It is noted that for 
a neutron target the % contribution is suppressed by the accompanying Dirac form factor 
pn (F™(t = 0) = 0), making this asymmetry sensitive to the CFF £ . However, in this 
case one expects a smaller single beam spin asymmetry that is also contaminated by other 
non-dominant CFF contributions. 
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Figure 11. KM 10 model predictions for DVCS single spin asymmetries: electron (4.13) [upper], 
transverse proton (4.14) [middle] and longitudinal proton (4.16) [lower] with E e xE p = 5 x 100 GeV 
(diamonds, squares) and E e x E p — 20 x 250 GeV 2 (triangles) EIC settings. Asymmetries versus 
azimuthal angle (j> for one selected bin (x B = 8.2 x 10" 3 , -t = 0.25 GeV 2 , Q 2 = 4.4 GeV 2 ) at 
5 x 100 GeV are shown in the left column for various GPD models: KM10 (dashed), KMlOa 
(dotted), and AFKM12 (solid). In the middle column the ^-dependence for the projection on the 
first non- vanishing harmonic is displayed for the same xb, Q 2 and beam energy values. In the 
right column the XB-dependencies is shown for the KM10 model prediction at —t = 0.25 GeV 2 , two 
different Q 2 values, Q 2 = 2.5 GeV 2 (filled squares and triangles, thick curves) and Q 2 = 13.9 GeV 2 
(empty squares and triangles, thin curves), and two different choices of beam energies, 5 x 100 GeV 
(red) and 20 x 250 GeV 2 (green). 



A single spin asymmetry measurement with a transversely polarized proton beam, 
cf. (2.33), provides another handle on the imaginary part of the helicity-flip CFF £. This 
asymmetry has in addition to the (j) dependence a <j) — (/)$ modulation. If the target spin in 
such a frame is perpendicular to the reaction plane (e.g., <j) — <j)s — 7r /2), the asymmetry 
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is dominated by a linear combination of % and £ CFFs. In the case that the target spin is 
aligned with the reaction plane (e.g., <j) — (j)$ — 0) the asymmetry 
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is formally dominated by a linear combination of CFFs % and £, cf. (2.20). In these 



asymmetries an additional relative kinematical factor J —t/AM^ appears. The middle row 
in Fig. 11 shows the sin(0 — (j)s) projection of the transverse proton beam spin asymmetry, 
which can also be rather large over a wide kinematical range. As in the case of the 
unpolarized cross section, discussed in the preceding section, this is caused by the fact that 
at smaller values of xb the "pomeron" behavior in % overtakes the kinematical suppression 
factors, see dashed and dotted curves. We may assume that such a "pomeron" behavior 
is also contained in £ . For our choice of ft sea = 1.5 the £ contribution will mostly cancel 
the % contribution, see (4.14) where F 2 (t = 0) « 1.79. In contrast to the electron beam 
spin asymmetry, for a neutron target the asymmetry is now more sensitive to the helicity 
conserving CFF %. For the cos(0 — <j)s) projection of the transverse proton beam spin 
asymmetry (4.15) the common expectation is that the parity-odd CFFs % and £ behave 
more gently at small xb and, hence, we expect that this observable is small in the EIC 
kinematics (not shown). 

Finally, we consider the longitudinally polarized proton beam spin asymmetry. Its 
projection on the dominant sin0 harmonic reads 



.sin^ V 1 -y I -t Ky o 



Fin-(fF 1 + ^F 2 )£ + f(F 1 + F 2 )n 



+ ••• . 

(4.16) 
It is sensitive to the imaginary part of CFF T~L and £, and other CFFs might contribute 
as well. As already noted, one expects that here the dominant CFF % behaves gently at 
small xb and models that incorporate such a behavior predict a rather tiny asymmetry 
(dotted and solid lines). In contrast, in the KM 10 model (dashed line), a rather big GPD 
H has been incorporated with a generic l/y / XB behavior at small xb- Hence, we get a 
sizable asymmetry for 5 x 100 GeV 2 beam energies which is getting smaller at higher beam 
energies 20 x 250 GeV 2 , see lower row on Fig. 11. We emphasize again that not much 
is known about the small-^B behavior of CFF %. We add that for a neutron target the 
asymmetry becomes sensitive to the CFF £ . 

Let us summarize the lessons from the approximated equations (4.13-4.16), quantified 
by numerics. The experimentally established 'pomeron' behavior of the CFF % predicts a 
large single beam spin and a large cos(0 — <j)s) projection of the transverse proton beam 
spin asymmetry for the EIC kinematics. If £ contains also a 'pomeron' behavior, the latter 
asymmetry can be weakened (amplified) for a positive (negative) imaginary part of £ . 
The remaining two single spin asymmetries cannot be predicted easily; however, based on 
common phenomenological/theoretical wisdom they are probably small. Let us note that 
the normalization of these asymmetries obviously depends also on the real part of the twist- 



35 



two associated CFFs and the remaining eight ones. As advocated above, a measurement 
of cross section differences are not affected by this normalization uncertainty. 

4.3 Further EIC opportunities 

An EIC machine provides further opportunities for DVCS studies: 

• Double spin flip experiments provide a handle on the real part of CFFs, however, in 
such measurements the spin-dependent BH cross section contributes. 

• As demonstrated by the HERMES collaboration, having a positron beam at hand 
allows also to separate the interference and DVCS harmonics in single spin target 
experiments. Measuring spin-dependent cross sections in the charge odd sector (in- 
terference term) and the charge even sector allows to extract CFFs from experimental 
measurements, based on minimal assumptions. 

• The large kinematical coverage of the proposed high-luminosity EIC (see Fig. 5) and 
the partial overlap with JLAB 12GeV kinematics raises the question: Can one utilize 
evolution, even at moderate xb values, to access GPDs away from their cross-over 
line? 

• Photon electroproduction off the neutron offers the possibility for a flavor separation. 

• Photon electroproduction off nuclei is a mostly unexplored experimental field. 

Below we will discuss a minimalistic version of the second point in more detail, namely, 
having an ^npolarized positron beam at hand. Let us mention here that a study of GPD 
evolution was presented in Ref. [119], however, we may conclude here that a wide coverage 
in Q 2 is extremely helpful in getting constraints on GPDs away from the cross-over line, 
however, a "measurement" of the GPD in the outer region certainly cannot be reached. 
DVCS on a "neutron target" is certainly needed for a GPD flavor decomposition. However, 
this program is more complicated than in DIS, since in the interference term the various 
CFFs are accompanied by nucleon form factors, see short discussions in the previous section. 
We will not discuss DVCS off nuclei, which is interesting in itself. It has been worked out 
theoretically for a spin-zero target, where one can adopt the equation from [136], and to 
some extent also for spin-one target [137-139], while the formalism for spin-1/2 nuclei can 
be adopted from the proton. 

We should also emphasize the EIC opportunities for Compton scattering measurements 
below the deeply virtual regime. 
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• Quasi-real Compton scattering can be measured over a rather wide energy range in 
anti-tagged electron scattering experiments, where the VCS cross section is peaked 
at Q 2 - 0. 

• We expect that at stage I binning of low photon virtualities, i.e., Q 2 < 1 GeV 2 , will 
be possible. 

Such measurements will provide understanding on the transition from the deeply virtual 
to the quasi-real regime. This, in turn, is needed if radiative electromagnetic corrections 
to photon electroproduction are to be elaborated in a more complete manner than they 
presently are. 

Finally, we should remind that other exclusive channels can be measured at EIC: 

• Deeply virtual production of light vector mesons can be employed for a partial flavor 
separation of quark GPDs. 

• J/ty production gives naturally access to the gluon GPD. 

• Also, experimental studies on deeply virtual production of pseudo scalar mesons, the 
production of two final meson states, time-like DVCS, and double DVCS may turn 
out to be feasible. 

We would like to add that deeply virtual production of light vector mesons and DVCS 
measurements at HERA collider experiments can be simultaneously described with a GPD 
formalisms [102, 123]. Whether the measurements, listed in the last item above, are actually 
feasible at EIC, can only be stated in terms of models. Thereby, based on phenomenological 
knowledge of the dominant GPD H, cross sections for time-like [141] and/or double [142- 
144] DVCS might be more or less realistically estimated, however, were not part of our 
studies. 

4.3.1 Uses of an unpolarized positron beam 

The isolation of the interference term, which contains the most valuable information on 
CFFs, is most easily done by forming charge asymmetries, which require a positron beam. 
We emphasize again, that the alternative Rosenbluth separation is expected to be more 
intricate and has not been so far either considered theoretically or explored experimentally 
(e.g., by the use of approximated expressions). Forming differences and sums of spin- 
dependent cross section measurements with both kinds of lepton beams allows to extract 
the pure interference and DVCS squared terms and might allow to quantify twist-three 
and gluon transversity effects. From such experiments one can extract the imaginary 
part of CFFs. Note that only an unpolarized positron beam is needed to perform such 
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a program for the single proton spin asymmetries - of course, for the projection of the 
single electron spin asymmetry a polarized positron beam would be needed. In double 
spin flip measurements one can use the same procedure to access the real part of the 
CFFs. Although existing data indicate that twist-three effects are small, as it is expected 
based on kinematic factors, the twist-three related CFFs are not necessarily small. Surely, 
one needs very high precision data to extract non-dominant twist-two CFFs or twist-three 
related ones. However, even obtaining only an upper limit is important for a determination 
of the systematic uncertainties of the (dominant) twist-two CFFs. 

Let us consider here only the lepton beam charge asymmetry (2.34) for an unpolarized 
proton. Its first harmonic is dominated by the real part of the twist-two related CFFs % 
and E , rather analogous to equation (4.13) for the electron beam spin asymmetry, 
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(4.17) 
It is shown in Fig. 12 that the GPD models predict a rather sizable lepton beam charge 




4> [rad] 



Figure 12. KM10 (dashed), KMlOa (dotted), and AFKM12 (solid) model predictions for the 
DVCS lepton beam charge asymmetry (2.34, 4.17) with E e x E p = 5 x 100 GeV 2 (upper row) and 
E e xE p = 20x250 GeV 2 (lower row). Left column: A c versus <\> for x B = 5.1xl0~ 3 , Q 2 = 4.4 GeV 2 , 
and t = -0.25 GeV 2 (upper panel) and x B = 5.1 x 10~ 4 , Q 2 = 4.4 GeV 2 , and t = -0.25 GeV 2 
(lower panel). Middle column: j^ s{ < n ^ amplitudes versus —t at same x B and Q 2 values for n = 1 
(filled diamonds and circles, thick curves), n = (empty diamonds and circles, thin curves), and 
n = 2 (stars). Right column: Dominant amplitude A^ s ^ ' versus x B for Q 2 = 2.5 GeV (thick 
curves) and Q 2 = 13.9 GeV 2 (thin curves) at t = -0.25 GeV 2 . 

asymmetries for both 5 x 100 GeV 2 (upper row) and 20 x 250 GeV 2 (lower row) beam 
energies. As for the electron beam spin asymmetry (upper row on Fig. 11), the predictions 
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from the KM 10 (dashed curves) and KM 10a (dotted curves) models are almost the same, 
illustrating that the CFF % contribution is rather unimportant in this observable (£ drops 
exactly out here), while the small deviation of AFKM12 model (solid line) indicate some 
sensitivity to CFF £. The sign of this asymmetry is governed by the effective 'pomeron' 
trajectory a p > 1. It has in the transition from the valence to the sea quark region a node 
(see upper right panel) with its position depending on —t. In the middle column we show 
besides the projection on the first even harmonic (filled diamonds and triangles) also the 
projection on the zeroth harmonic (empty diamonds and triangles), which is dominated by 
twist-two associated CFFs and, thus, (ant i) correlated with the first harmonic. The second 
harmonic (stars) is sensitive to twist-three associated CFFs, which are set here to zero, and 
also depend on the twist-two associated CFFs. Therefore, the latter induce only a small 
deviation, compatible with zero within one standard deviation. 

5 Partonic interpretation at small xb 

Intensive GPD studies (up to NNLO accuracy) of small-^B DVCS data measured by HI and 
ZEUS collaborations have been performed, where it turned out that the functional form 
of the ^-dependence cannot be pinned down and an access to the CFF £ is not feasible 
when having only unpolarized DVCS cross section and the lepton beam charge asymmetry 
measurements [6] available [110]. A high-luminosity EIC experiment with transversely 
polarized protons certainly provides the opportunity for precise measurements of CFFs 
and to explore their partonic interpretation in the small-xe region, i.e., xb < 0.01. As 
we argued in Sect. 4, the set of relevant twist-two associated CFFs is then reduced to H 
and £ only, and, moreover, valence quark contributions can be safely neglected. From 
our discussion there it is also obvious that the former assumption, which is used now, 
can be experimentally cross-checked. We will also use the fact that the real parts of the 
remaining two CFFs is locally tied to their imaginary parts, see (4.1), which is implemented 
in our GPD model and is in fact a more general consequence of the dispersion relation and 
the effective 'pomeron' behavior. Hence, we can restrict ourselves to two observables, 
namely, the unpolarized DVCS cross section (4.11) and the single transverse proton beam 
asymmetry (4.14), which now simplify to 
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In the partonic interpretation of DVCS data we are in the first place interested in the 
transverse distribution of sea quarks and gluons at small xb for an unpolarized and for a 
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transversely polarized proton. In Sec. 5.1 we explore by least-squares fitting the extrac- 
tion of both GPD H and E from the aforementioned observables at stage II of an EIC. 
In Sect. 5.2 we present a detailed study of the extraction of transverse polarized parton 
distributions. We also perform there the Fourier transform of our GPD model fit results 
to the impact space, where experimental uncertainties are propagated and extrapolation 
errors are taken into account. Finally, in Sect. 5.3 we discuss the importance of such a 
measurement for the qualitative understanding of the proton spin decomposition. 

5.1 Extraction of GPDs H and E from high energy EIC pseudo data 
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Figure 13. A model dependent extraction of GPD H from cross section measurements of the 
HI /ZEUS collaborations (left) and from a combined fit that includes EIC pseudo data (right) with 
beam energies E e x E p = 20 x 250 GeV 2 . The HERA collider data are taken from [5] (circle), [3] 
(squares, diamonds), and [6] (triangle- up, triangle-down, rectangle). 

To illuminate how GPDs are experimentally constrained in the small- x region at 
present, we will present here also new fits to the world DVCS data set at large W (small 
xb) that includes the propagation of experimental uncertainties. At the HERA collider 
experiments HI [1, 3, 4, 6] and ZEUS [2, 5] the unpolarized DVCS cross section could be 
measured at large W, — t < 1 GeV 2 , and with a large lever arm in Q 2 . Although ~ 200 
data points were published, we only consider 35 as statistically independent. Apart from 
those for the differential cross section that are displayed on the left panel on Fig. 13 versus 
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N 


a 


a' [GeV" 2 ] 


b [GeV -2 ] 


S2 


s A 


Po a 


0.152 


1.158 


0.100 


2.800 


0.513 


-0.210 


,^sea 


- 


- 


0.090 


2.858 


0.508 


-0.208 


Sp sea 


- 


- 


0.009 


0.035 


0.038 


0.011 


Po 


(0.448) 


1.247 


0.100 


2.000 


-4.806 


1.864 


p G 


- 


- 


0.063 


2.086 


-4.739 


1.835 


5p G 
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0.088 


0.163 


0.212 


0.106 


E 


hi 


a 


a' [GeV" 2 ] 


b [GeV~ 2 ] 


S2 


sa 


^sea 


1.500 


1.158 


0.020 


2.800 


0.513 


-0.210 


^sea 


1.451 


1.164 


0.023 


2.779 


0.524 


-0.213 


Sp sea 


0.307 


0.005 


0.012 


0.049 


0.104 


0.026 


Po 


(-0.51) 


1.247 


0.050 


2.000 


-4.806 


1.864 


p G 


(-0.49) 


1.295 


0.001 


1.961 


-4.687 


1.803 


Sp G 


(0.06) 


0.216 


0.252 


1.092 


0.048 


0.077 



Table 2. AFKM12 model parameters (po) and their fitted values (p) together with standard 
uncertainties (Sp) for sea quarks (superscript sea ) and gluon (superscript G ) components of GPDs 
H and E at the input scale Q 2 , = 4 GeV 2 . The values in parentheses are fixed by sum rules. 

—t we also included the following ^-integrated cross section measurements from ZEUS 

7.5 GeV 2 < (Q 2 > < 85 GeV 2 (W) = 89 GeV [2], (5.3) 

7.5GeV 2 < (Q 2 > < 70GeV 2 (W) = 104GeV [5]. 

In our LO fits to these data we can only ask for the sea quark and gluonic components of 
GPD H, where both components can be separated to some extent due to the large Q 2 lever 
arm [110]. Since the experimental uncertainties are large, the functional form of the —t 
dependence for sea quark (and gluon) GPD H cannot be determined by x 2 /d.o.f . w 1 model 
fits, done here at LO with a dipole (dashed) or an exponential (dash-dotted) residual t- 
dependence. Thereby, also the "pomeron" slope parameter a f for the sea quark (and gluon) 
content cannot be determined. 

To explore the potential of the EIC measurements at stage II, we use in the following 
pseudo data for the unpolarized DVCS cross section (5.1) and the transverse target spin 
asymmetry (5.2) for the beam energies 20 x 250 GeV 2 , as specified in Sect. 3. Thereby, 
we utilized the flexible AFKM12 model, introduced in Sect. 4, in which GPD H and 
E have a different "pomeron" slope parameters, see po values in Tab. 2, but the same 
residual ^-dependencies. This choice guarantees that positivity conditions for GPDs at 
zero skewness are mostly satisfied [145]. The experimental uncertainties were estimated 
as before (statistical uncertainties from the MILOU simulation, which are rescaled for the 
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DVCS cross section, 5% systematic uncertainty on cross section level, 3% uncertainty of the 
BH cross section in the subtraction procedure (4.8), and 5% beam polarization uncertainty). 
The exponential ^-dependence of the CFFs drastically increases the subtraction uncertainty 
at large — £, and the net uncertainty in this region can become very large at larger y values 
(lower xb values in particular at low Q 2 ), see right panel on Fig. 13. Since for the GPD E 
we took a model with a positive ft sea = 1.5, the transverse target asymmetry is becoming 
small. Fig. 14 shows the pseudo data for this asymmetry together with the model curve 
(solid), used to generate the asymmetry. Also shown is the prediction from an analogous 
model which, however, has a negative ft sea = —1.5 value (dashed curves) and one with 
vanishing CFF £ (dash-dotted curves). Certainly, the predictions of all these three models 
are experimentally distinguishable. 
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Figure 14. EIC pseudo data (diamonds) for the transverse target spin asymmetry (5.2) at beam 
energies E e x E p = 20 x 250 GeV 2 are shown together with AFKM12 GPD model predictions, where 
GPD E sea is taken as large positive (solid), vanishing (dot-dashed), and large negative (dashed), 
respectively. 

We performed a combined least-squares fit to the EIC pseudo data for the unpolarized 
DVCS cross section and the single transverse proton beam asymmetry together with the 
HERA collider measurements, shown in Fig. 13 and (5.3). Altogether we included 2732 
data points, where the EIC pseudo data were generated as specified above and contain 509 
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data points for the unpolarized cross section, obtained from 21 {Q 2 , xsj-bins: 

3.16GeV 2 <Q 2 < 5.62 GeV 2 , 2.5 x 10" 4 < x B < 1.0 x 10" 2 (8 bins), 
5.62 GeV 2 < Q 2 < 10.00 GeV 2 , 4.0 x 10~ 4 < x B < 1.0 x 10~ 2 (7 bins), (5.4) 
10.00 GeV 2 < Q 2 < 17.78 GeV 2 , 6.3 x 10" 4 < x B < 1.0 x 10" 2 (6 bins), 

and 2188 data points for the single transverse beam spin asymmetry as function of 0, 
obtained from 24 {Q 2 , x B }-bms: 

3.16 GeV 2 <Q 2 < 5.62 GeV 2 , 1.58 x 10" 4 < x B < 1.0 x 10" 2 (9 bins), 
5.62 GeV 2 < Q 2 < 10.00 GeV 2 , 2.51 x 10~ 4 < x B < 1.0 x 10~ 2 (8 bins), (5.5) 
10.00 GeV 2 <Q 2 < 17.78 GeV 2 , 3.98 x 10" 4 < x B < 1.0 x 10" 2 (7 bins). 

In the fit we released all 19 model parameters, which are partially correlated. In particular, 
the normalization factor ft sea of GPD E is strongly correlated to the skewness parameters S{. 
Obviously, the hypothesis of a dipole ^-dependence yields an unacceptably large x 2 /d.o.f. 
value, while the exponential ansatz provided, as it should, almost the textbook value of 
one, x 2 /d.o.f . = 0.97. The extracted parameters and their standard uncertainties are listed 
in Tab. 2. The slope parameter of H sea can be well extracted with less than two standard 
deviations away from the input model parameter value. The normalization of this GPD for 
fixed PDF parameters is also rather robust, as indicated by small deviations of extracted 
skewness parameters from the model parameters. Since the pseudo data constrain the 
^-dependence, the correlation of normalization parameters and t-slope (or dipole mass) 
parameters is much less pronounced than in the fits to the HERA collider data. For the 
GPD H G the uncertainty for ol is of the order of its model parameter value 0.1 and the 
relative uncertainty of the residual ^-dependence is now of the order of 7% rather than 1% 
as for sea quarks. The relative size of skewness parameter uncertainties for gluons is on 
the same 5% level as for quarks. For GPD £" sea the "pomeron" intercept, normalization 
K sea and skewness parameters are well reproduced by the fit, where the ft sea uncertainty is 
of the order of 20%. The moderate size of this uncertainty also reflects the correlation of 
the normalization with the skewness parameters, where the latter is now more than twice 
larger than for GPD i7 sea . The uncertainties for the t-slope parameters are only about 40% 
larger than for GPD i7 sea and are still reasonably small. For E G already the "pomeron" 
intercept parameter has a very large uncertainty, which will induce a huge normalization 
uncertainty. Note also the t-slope parameters have big uncertainties, and they are also 
correlated with the remaining ones. In general we found that with our conservative fitting 
strategy it is impossible to access the gluonic component of GPD E from the employed set 
of DVCS pseudo data. It is a standard procedure to reduce the set of parameters to those 
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that are not strongly correlated. This will also reduce the size of uncertainties, however, 
certainly one should bear in mind that this procedure increases the theoretical bias. 
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Figure 15. Least-squares fit extraction of sea quark GPD H sea (left) and gluon GPD H G (middle) 
from a dipole ansatz (gray area surrounded by dashed curves) and an exponential ansatz (gray area 
surrounded by dashed-dotted curves) using only the HERA collider data. The results of a combined 
HERA/EIC fit including pseudo data for the unpolarized DVCS cross section, c.f. Fig. 13, and the 
transverse target spin asymmetry A^x , c -f- Fig. 14, using an exponential ansatz are shown as 

light orange area (surrounded by solid curves). In addition for the first time the sea quark GPD 
^sea cou y j^ extracted (right panel). 

In Fig. 15 we compare the resulting GPDs from fits to the HERA data alone and to 
the combined HERA+EIC data at Q 2 = 4 GeV 2 , xb = 10 -3 , and variable — t (covering 
the HERA region). In the right panel one realizes that the uncertainty of the sea quark 
GPD i7 sea 5 which is to certain extent constrained by HERA data, can be strongly improved 
in particular at smaller —t values. The gluon GPD H , displayed in the middle panel, 
is extracted by means of the Q 2 evolution and it is rather weakly constrained by HERA 
DVCS data only. Here the inclusion of stage II EIC data yields a large improvement, even 
if the used lever arm in Q 2 , compared to HERA kinematics, is still rather limited. As 
emphasized above, information on the GPD E can only be obtained from a new lepton- 
proton scattering experiment with a transversely polarized proton beam. In the right panel 
it is clearly demonstrated that the sea quark component of this GPD can be extracted with 
relatively small uncertainties. As explained above, from the utilized pseudo DVCS data the 
gluonic component of GPD E (not shown) cannot be reliably accessed using our flexible 
GPD models. 

5.2 Transverse spatial imaging 

One of the main goals of GPD phenomenology is to provide the transverse spatial distribu- 
tions of partons as function of the momentum fraction x. The simplest proposal to obtain a 
rough idea of such parton distributions is based on the Fourier transform of the amplitude 
[48]. Following the common experimental procedure one would extract the ^-dependence 
from a fit to a given (xb, Q 2 ) bin. Utilizing the HERA data for DVCS and exclusive J/ip 
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production and saying that the former process is quark dominated while the latter is gluon 
dominated, one immediately concludes from the experimental findings that the exponential 
t-slope parameter for DVCS cross section is larger than for the J/ifj cross section, 

Sdvcs « 6 GeV" 2 > Bj^ « 4 GeV" 2 , 

meaning that sea quarks are more spread out in transverse space than gluons. However, 
we may note that this rather generic interpretation is based on the assumption that the 
proton helicity non-conserved CFFs and/or amplitudes play no important role and that 
skewness effects are unimportant. We also emphasize that in a partonic interpretation the 
accessible lever arm in —t is restricted by the DVCS requirement — t <C Q 2 , which ensure 
that possible higher twist contributions, twist-four and higher, are small. 

To quantify possible differences between a GPD interpretation and the aforementioned 
procedure, our GPD fit result from the preceding section is compared with (half of) the 
exponential t-slope of the differential DVCS cross section. The latter is extracted by fits 
to the ^-dependence in a given {xb, Q 2 }-bin of the pseudo data (5.4) by means of the 
exponential model 

1 dt W =n,exp{2M}, (5.6) 

where ni and hi are the two fitting parameters used in bin i. The x 2 /d.o.t value in these 
fits is usually around one, where the propagated standard error can be rather large for 
the lowest XB-bins, due to the BH subtraction procedure. From our GPD fit, we employ 
both the (sea) quark GPD H on the cross-over line and the square root of the predicted 
differential DVCS cross section, the latter containing also additional contribution due to the 
non- vanishing CFF £. ^From both of these quantities we calculate an effective exponential 
^-slopes in the same manner. E.g., for the GPD on the cross-over line such a t-slope reads 

6 eff (x B , Q 2 ) = -±- In ^f%\ (5-7) 

h - h H(x,x,ti, Q z ) 

with ti = -0.03 GeV 2 , t 2 = -1.5 GeV 2 , and x = x B /{2 - x B ). 

In the left panel of Fig. 16 we show the results of the local exponential model (5.6) fits to 
the DVCS cross sections versus x B for Q 2 = 4.1 GeV 2 , Q 2 = 7.3 GeV 2 , and Q 2 = 12.9 GeV 2 
as empty triangles-down, circles, and triangles-up, respectively. The effective t-slopes (5.7) 
of the GPD on the cross-over line, extracted from the combined GPD fit, are shown as 
thick solid, dash-dotted, and dashed curves, respectively, while effective t-slopes of GPD- 
flt-predicted cross sections are plotted as thin lines. Note that thick curves contain the 
uncertainty bands, arising from the propagation of pseudo data uncertainties. As one 
realizes, in our large E scenario the t-slope of the DVCS cross section (5.6) is relatively flat 
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Figure 16. Exponential t-slope parameters from the model fit (5.6) to the DVCS cross section 
(empty symbols) and of ^siaH from the model fit (5.8) to cross section and asymmetry data (filled 
symbols) as well as effective t-slope parameters (5.7), obtained from our GPD fit, for the GPD H 
on the cross-over line (thick lines), the DVCS cross section prediction (thin lines), and the zero- 
skewness GPD H (error bands in the right panel) are displayed as function of xb for three different 
Q 2 values: 4.1 GeV (triangles-down, solid curves, blue band), 7.3 GeV (circles, dash-dotted curves, 
green band), and 12.9 GeV 2 (triangles- up, dashed curves, red band). Uncertainty bars arise from 
local fits to EIC pseudo data, while curves and uncertainty bands originate from a combined GPD 
model fit to HERA collider and EIC pseudo data. 

w.r.t. both xb- and Q 2 -dependence, see empty symbols. This behavior differs drastically 
from that of the extracted GPD H (thick curves), which has in addition also larger slope 
values. Nevertheless, the 6-slope values extracted via exponential fit (5.6) are consistent 
with the effective DVCS cross section t-slope (thin curves), obtained from our GPD fit, and 
evaluated analogously to (5.7). In fact, the DVCS cross section (5.1) is in our model given 
as linear combination of two exponentials with different slope parameters. The behavior 
of GPD H (left panel: thick curves) is partially compensated by the appearance of CFF 
£ . Note that for a small/vanishing GPD E scenario these differences in the 6-slope values 
would die out and the differences in the ^-dependence of the CFF modulus \H\ and QrnH 
can be considered to be small in the studied kinematical region. Consequently, under these 
circumstances and restricting the fits to LO accuracy, one can extract the ^-dependence of 
GPD H(x,x,t, Q 2 ) directly from the DVCS cross section measurements. 

To go beyond the % dominance hypothesis in such cross section fits, one can utilize 
measurements of the single transverse proton spin asymmetry (5.2). We recall that this 
asymmetry is sensitive to CFF £ and so the appropriate strategy is to use cross section and 
asymmetry data simultaneously in an analysis. To perform local fits for fixed xb and Q 2 
one may set the real part of CFFs to zero. Alternatively, one can utilize a Regge-inspired 
ansatz, e.g., rather analogous as in (4.1), and perform fits for given Q 2 , see, e.g., Ref. [146]. 
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We performed such 2x4 parameter (n, a(0), a! , 6 for CFFs % and £) fits in the three Q 2 
bins of the pseudo data sets (5.4) and (5.5). To compare the propagation of uncertainties 
with our local DVCS cross section fits, we used the extracted values for the two Regge 
trajectories a + a't in local 2x2 parameter fits to the EIC pseudo data. The 6-slope 
parameters of the imaginary parts are defined as 

$smHi(t) — rti exp{b{t} and 9m^(t) — rii expj&it} , (5.8) 

and the slightly ^-dependent phase for a given {xb, Q 2 }-bin i was considered to be known. 
Furthermore, we restricted the set of asymmetry data (5.5) to those of the DVCS cross 
section (5.4). The results for the exponential t-slope parameter of SsmH are presented on 
the left panel in Fig. 16 as filled symbols for Q 2 = 4.1 GeV 2 (triangle-down), Q 2 = 7.3 GeV 2 
(circle), and Q 2 = 12.9 GeV 2 (triangle-up) and, as expected, they are compatible with the 
effective t-slope of the GPD H{x 1 x 1 t 1 Q 2 ) 1 extracted from our GPD fit (thick curves). 
However, the propagated uncertainties in these local fits are larger than in the previous 
ones, reflecting the fact that, particularly at larger xb, the asymmetry uncertainty can 
get large, see Fig. 14. Surely, assuming that both the assumed uncertainty distribution 
(Gaussian) and the model is correct, the uncertainty propagation in global fits, e.g., with 
a Regge-inspired ansatz, provides a much smaller error. It should be noted that both of 
these assumptions are only true with a certain probability. 

Next we consider the effective t-slope, analogously defined as in (5.7), of the unpolarized 
quark GPD without skewness dependence, 

q(x,t,Q 2 ) = H(x,7 1 = 0,t,Q 2 ), (5.9) 

which is of great interest with regard to the transverse distribution of quarks. The result 
of our model is presented on the right panel of Fig. 16 and it can be compared to slopes of 
the GPD H(x, rj = x, t, Q 2 ) on the cross-over line in the left panel (thick lines). First it is 
observed that with growing Q 2 not only the x-slope 10 of the effective t-slope decreases as 
it is also seen for GPD H(x,x,t, Q 2 ), but also its intercept drops. Loosely spoken, such a 
behavior means that both the "pomeron" slope parameter a! and the value of the residual 
t-slope parameter decrease. This behavior of the zero-skewness GPD looks more natural 
and it is naively expected from the double log asymptotic behavior also for the GPD on the 
cross-over line [147]. However, with our AFKM12 model we illustrate by the thick curves 
in the left panel of Fig. 16 that specific choices of skewness parameters combined with t- 
slope parameters can provide also a rather flat effective residue dependence and a stronger 



10 Since x ~ xb/2 is valid at small xb, twice of the XB-slope that can be read off from Fig. 16 can be 
equivalent ly considered as x-slope. 
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decrease of the "pomeron" slope parameter a! . Since the 77 —¥ limit commutes with 
the Q 2 -evofution, we conclude that evolution entangles the ^-dependence with skewness 
dependence. Hence, the factorization of the t and skewness dependence, assumed in our 
model at the initial scale, does not hold true under evolution (otherwise the effective 
slope parameters should evolve similarly). The reader might be also surprised that the 
uncertainty bands of the effective slope for the GPD on the cross over line remain tiny 
(thick curves on the left panel), while those in the forward (zero-skewness) case get sizable 
with increasing Q 2 and decreasing x. This is caused by a naive truncation of the covariance 
matrix, i.e., removing rows and columns belonging to the skewness parameters 52 and 54 
in the forward limit, which also alters the (anti)correlation of uncertainties that ensure 
the smallness of the uncertainties for GPD H{x 1 x 1 t 1 Q 2 ). It is beyond the scope of this 
paper to study the uncertainties that arise from the extrapolation to 77 — >► in more depth, 
however, once high precision data will become available, one should also worry about the 
model bias in the extrapolation of GPD F(x, 77, t, Q 2 ) from 77 = x to 77 = 0. 

Finally, we would like to illustrate that the EIC measurements of the ^-dependence 
in the region — t < 1.5 GeV can provide a probabilistic interpretation of the transverse 
distribution of sea quarks and partially also for gluons. The Fourier transform of the 
zero-skewness GPD (5.9) into the impact parameter space, 



q(x,b,Q 2 



JfJ^e- iKrb H(x, V = 0,t = -K 2 ,Q 2 ) (5.10) 

= — d\t\J [b^\j H(x, rj = 0, t, Q 2 ) , 

provides in the infinite momentum frame the probability of scattering on a quark as a 
function of its momentum fraction and transverse distance b —\b\ from the proton center, 
where 1/Q is considered as the resolution scale [68]. Since for an unpolarized struck quark 
and proton, no direction in the transverse plane is preferred, the integration over the polar 
angle in (5.10) yields a Bessel transform (Jk(x) denotes the Bessel function of order k) 
and results in a parton density is symmetric under rotation of the two-dimensional impact 
parameter vector b. For a transversely polarized proton, e.g., the polarization vector is 
pointing in the x direction, one finds that the parton density is given by the unpolarized 
one (5.10) and a distortion in y direction that is governed by the strength of GPD E, 

q\x, b, Q 2 ) = q(x, b, Q 2 ) - J^^Le(x, b, Q 2 ) (5.11) 



(x,ri = 0,t,Q 2 



= — / d\t\ Jo(bJ\t\) H + ^-^JAbJki) E 
4vr,/o V v ' 7 2bM p V v ' 7 

Before we present the resulting parton densities (5.10) from the combined GPD model 
fit to HERA and EIC pseudo data, let us shortly discuss the peculiarities in the uncertainty 
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estimation. The uncertainty of the resulting parton densities is, besides the propagated 
experimental uncertainties, also dictated by the possible uncertainties caused by extrap- 
olations from the accessible kinematical region, namely, (i) extrapolation of the skewness 
parameter dependence 77 = x to 77 = 0, discussed above, (ii) extrapolation of ^-dependence 
from the experimental minimal — t value —t\ to —t = 0, as well as (iii) from maximal acces- 
sible value —£2 to — t = 00. These rather intricate extrapolations are fortunately governed 
by the boundary condition, 

r poo 

q(x,Q 2 )=H(x,r] = 0,t = 0,Q 2 ) = d 2 b q(x,b,Q 2 ), (5-12) 

J J —CO 

arising from the reduction of GPD H in the kinematical forward limit to the standard 
unpolarized PDF q. Hence, the normalization of the (integrated) parton density (5.10) 
is also entirely determined by the PDF normalization. To simplify our study, we restrict 
ourselves to Q 2 = 4GeV 2 , where in our model the t- and skewness dependencies factorize, 
as discussed above and exemplified also by the agreement of the effective slope parameters 
in the r\ — x and 77 = case, see thick solid curves on the left and right panels on Fig. 16. 

A model analysis studying the challenges of extrapolation in —t beyond the exper- 
imentally accessible range has been presented for the differential cross section in [148] 
and we essentially agree with the conclusion that with an EIC imaging is feasible for 
0.1 fm < b < 1.5 fm (or even in a wider range). Let us add some mathematical insight and 
let us point out methods to increase the quality of the extrapolations. With our model 
hypothesis the ^-dependence of the zero-skewness GPD is essentially constrained by the 
EIC pseudo data in the region 0.03 GeV 2 < -t < 1.5 GeV 2 . 

The uncertainty of the extrapolation into the region [0, — 1\) is associated with the 
contribution 

A iq (x, 6, Q 2 ) = ±- f 'd\t\ J (&v1*[) H(x, V = 0, *, Q 2 ) , (5.13) 

from which one can easily obtain estimates. Although q(x : Q 2 ) = H(x, 77 = 0, t, Q 2 ) at t = 
is very well known, which makes this an interpolation problem rather then an extrapolation 
one, let us here calculate A±q from the knowledge of H in the vicinity of —t\ by making 
use of a truncated Taylor series where the uncertainty is equated with the remainder. 
Consequently, to first order accuracy we have 

A iq (xAQ 2 ) » ^ J i( b \VW\) H(x,0,t 1 ,Q 2 ) + ^^J 2 {b\^\h\) A-H(x,0,h,Q 2 " 

2|fl|3/2 7 (b fW~\\ d " 



S iq (x,b,Q 2 ) ~ ^Jal^vW) ^H(x,0,h,Q 2 ), (5.14) 
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where the derivative of H(x, 0, ti, Q 2 ) can be evaluated numerically. For small —t\ we can 
roughly estimate the value of the second order derivative in terms of the transverse width 

(b 2 )(x,Q 2 )=4j t \nH(x, V = 0,t,Q 2 )\ , (5.15) 



t=0 



e.g., for a p-pole ansatz we find 



^H(x,0,h, Q 2 ) « l -^(b 2 ) 2 {x, Q 2 )q(x, Q 2 ) , 

where the result for an exponentially functional form in t follows from the limit p — >► oc. 
For a realistic value of (b 2 ) ~ 0.35 fm in the small x region, we find that the extrapolation 
uncertainty is of the order of 10 -4 in units of q(x, Q 2 )/fm 2 . However, this uncertainty 
becomes important in the large b > 1/y^ftij ~ 1 f m region, dominated by the contributions 
from the small —t region. Model analyzes provide a relative uncertainty on permill level for 
6^1 fm, which, however, will increase to the few percent level for b « 1.5 fm and will then 
grow fast for increasing b. Hence, with our EIC pseudo data we can resolve the transverse 
distribution up to a distance of ~ 1.5 fm. To estimate the uncertainty of the extrapolation 
into the {—£2, °°] region we naively use 

1 r°° 

A 2 q(x, b, Q 2 ) = ~r d\t\ Jo(b^/\t\) H(x, Tj = 0,t, Q 2 ) . (5.16) 

It is most important for the small-6 region, where in particular for b — we have 

1 f°° 
A 2 q(x,b = 1 Q 2 ) = — d\t\H(x,7 1 = 0,t,Q 2 ). (5.17) 

47r J\t 2 \ 

The relative uncertainty at b — is easily evaluated and we find for an exponential or 
p-pole form 

A 2 q(x,b = 0, Q 2 ) = e(b 2 )t2/4 ^ A 2 g(x,6^0,Q 2 ) = / (6 2 )t 2 ^ pfl 



g (x,6 = 0,Q 2 ) q(x,b = 0, Q 2 ) V 4 P 

respectively. Assuming an exponential functional form it is with —t<i « 1.5 GeV 2 well 
under control and results in a ^ 3% correction, which, however, would increase for a 
dipole form to ~ 40%. To reach the 10% accuracy level, one must increase —t<i ~ 8 GeV 2 , 
which requires a big Q 2 value to ensure the validity of DVCS kinematics. Fortunately, the 
error, e.g., for 0.1 fm < 6, gets already on the 10% level for —t 2 ~ 3.5 GeV 2 . Under these 
circumstances, one may rely on extrapolation techniques, e.g., based on conformal mapping 
or Pade approximation, to minimize the uncertainty. Note also that the uncertainty of 
extrapolation into the {—£2,00] region may be also associated with a relative uncertainty 
that grows fast with increasing b. In the following the uncertainty is calculated according to 
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(5.17) and estimate numerically by assuming two alternative hypotheses, namely, that the 
^-dependence falls off exponentially or with 1/t 2 , where for a given b value always the larger 
uncertainty is taken. For simplicity we will neglect the uncertainty from the extrapolation 
(interpolation) into the region {— £i,0], which is entirely justified for b < lfm and as it 
would be hardly visible in the visualization of the parton densities for b < 1.5 fm. Finally, 
the uncertainty from the extrapolation into the large —t region was added in quadrature 
to the one propagated from the (pseudo) data. 
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Figure 17. Parton densities at x = 0.001 and Q 2 = 4GeV versus impact parameter b were 
obtained from a combined least-squares fit to the HERA collider and EIC pseudo data: relative 
densities (lower row) and their values at b x = for the unpolarized sea quark parton densities of 
a unpolarized proton (left), a transversely polarized proton (middle), and the unpolarized gluon 
parton density of a unpolarized proton (right), its value is rescaled by a factor 0.19. 

In the left and right columns on Fig. 17 the sea quark and gluon parton densities (5.10) 
at x — 10 -3 and Q 2 = 4 GeV 2 are shown as a relative density plot versus b y and b x (lower 
panels) and for b x — as function of b y (upper panels). Note, the gluon density is rescaled 
by a factor of 0.19. Since for sea quarks the propagated uncertainty is small and it was 
assumed that there is no cross-talk between t- and skewness dependencies, and the PDF 
uncertainties [fixed PDF parameters, see also boundary condition (5.12)] are neglected, the 
final uncertainties are smaller than for the corresponding ^-dependent GPD at the cross- 
over line shown in the left panel of Fig. 15. Nevertheless, the increase of the uncertainty 
in the vicinity of b — due to the {— £2, 00] extrapolation is visible for the quark density 
(upper left panel) and much more pronounced for the gluon density (upper right panel). 
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Generally, the larger error for the gluon density is mainly based on the fact that the DVCS 
process alone does not allow to pin down this quantity on the same quantitative level as 
for sea quarks, see earlier discussions and Fig. 15. We again emphasize that the functional 
form of the ^-dependence will influence the uncertainties related to the extrapolation error. 
For instance, a power-like falloff will increase the DVCS amplitude in the accessible large 
—t region and therefore decrease the experimental uncertainties in this region, however, on 
the other hand, the uncertainties of the extrapolation into the {—£2, °°] region will become 
more important. 

Apart from the uncertainties that appear in the unpolarized parton densities, we also 
have the normalization uncertainty of GPD E, which is not protected by a boundary 
condition. Nevertheless, we found in our model that this uncertainty is not large. This is 
illustrated in the middle column of Fig. 17, where we display the sea quark density for a 
transversely polarized proton. 

5.3 Angular momentum sum rule 

Finally, we shortly discuss the role of the EIC measurements in elucidation of the Ji spin 
sum rule [71]. This rule states that the proton spin 

1= Y, J q (Q 2 ) + J G (Q 2 ) (5-18) 

q=u,d,s,--- 

is built from quark and gluon angular momenta J q and J , which are defined via a gauge 
invariant decomposition of the QCD energy momentum tensor. Note that several other 
decompositions have been proposed, which are related to Ji's ones by reshuffling a certain 
amount of angular momentum fraction A J, i.e., 

J2<J q ^J2 jq + AJ and J°^J G -^J, (5-19) 

q q 

where one may take the freedom to define A J as the expectation value of a gauge variant 
operator in order to reach a partonic interpretation of the gluon component in terms of spin 
and orbital angular momentum, e.g., to arrive at the Jaffe-Manohar spin sum rule [149]. 
Ji's decomposition implies that the partonic components of the proton spin are given by 
the momentum fraction part, called here A, and the anomalous gravitomagnetic moment 
B, 

J* (Q 2 ) = l -A\Q 2 ) + ^(Q 2 ), J A B J' (Q 2 ) = J\x x{*}\x,ri = 0,t = 0, Q 2 ) , 

(5.20) 
which are given by the first moments of GPDs H and E, respectively. A phenomenological 
quantification of this sum rule is a highly intricate task, which is often trivialized by entirely 
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relying on simple-minded GPD (or even transverse momentum dependent PDF) models or 
assumptions. 

The definitions in Ji's angular momentum sum rule allow to employ any other QCD, 
i.e., field theory based framework to quantify the quark and gluon angular momenta. Most 
promising for achieving this goal are lattice gauge field simulations and once reliable results 
can be obtained, for a review see [150], one may incorporate them in GPD models. There 
are various systematic uncertainties in the lattice estimation of angular momentum carried 
by sea quarks and gluons and hence their phenomenological determination is an important 
task for the future. 

As we have seen, DVCS measurements at an EIC will allow to access the GPD E at the 
cross-over line and allows, in a model dependent manner, to extract also its normalization 
in the forward kinematics. In fact, what we called anomalous magnetic moment of sea 
quarks is simply related to their angular momentum: 

J sea = -(l + /^ sea )A sea , 

where the phenomenological value of the momentum fraction is at Q 2 = 4 GeV given by 

A sea (Q 2 = 4GeV 2 ) wO.15. 

6 Summary 

We show on some selected examples the physics case for DVCS measurements at a proposed 
EIC. Pseudo data were generated by the MC program MILOU that is tuned to HERA 
collider DVCS measurements. Full detector simulations have not yet been included; it 
was rather assumed that the systematical uncertainty for cross section measurements is 
on the 5% level. The statistical uncertainties of these simulations have been included in 
model predictions for various single spin asymmetries, electron charge asymmetries, and 
unpolarized cross sections, covering the EIC kinematics at stage I and II. 

We illustrated that present GPD models, constrained by global fits to present DVCS 
data, provide a variety of EIC predictions, where in particular the t-dependence of the 
different models is poorly known and can be constrained to a large degree by EIC measure- 
ments. We did not discuss in completeness the extraction of CFFs, which can be done by 
having a polarized positron beam at hand. Here already a unpolarized one would help to 
have a cleaner access to twist-two associated CFFs. In particular, it can be used to isolate 
the interference term, which contains the most valuable information. This also provides 
an experimental cross-check for the smallness of the ^-integrated interference term in cross 
section measurements. Rosenbluth separation, as it was worked out in Sect. 2.1, provides 
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another handle on the isolation of BH and DVCS cross sections. This technique has to be 
explored further for the access of higher harmonics in the interference term. At present it 
is not known to what extent this method can be employed in a model independent manner, 
however, certainly it looks more intricate than in the case of unpolarized DIS or elastic 
form factor measurements. 

While the access to CFFs and GPDs at lower beam energies requires the measurement 
of many observables the situation becomes simpler at higher energies. Here, we can assume 
that only two twist-two associated CFFs % and £ show a "pomeron" behavior and are as 
such accessible in these kinematics. Moreover, their real parts are small compared to the 
imaginary and their phases are approximately given by an effective "pomeron" trajectory. 
Therefore, they can be accessed by a measurement of the DVCS cross section and the 
single transverse proton spin asymmetry. Thereby, at large electron energy loss y and —t 
the DVCS cross section may drop drastically and perhaps cannot be measured. In such 
circumstances one can use the single electron beam spin asymmetry measurements, which 
are predicted to be sizable at large y. 

For parton imaging, it was illustrated that in a large GPD E scenario its ^-dependence 
extracted from the unpolarized cross section does not necessarily match the ^-dependence 
of the sea quark GPD H sea . To extract in such a situation the ^-dependence of CFFs % 
and £ one may neglect the real part or use global fits, e.g., with a Regge-inspired ansatz 
in a given Q 2 bin. Taking EIC pseudo data generated from MILOU and propagated to 
the predictions of AFKM12 model, we studied the error propagation to the sea quark 
and gluonic components of GPDs H and E by means of least-squares fits. Thereby, it 
turned out that with our rather flexible model the sea quark component of both GPDs can 
be pinned down quite precisely, while the knowledge of the gluon GPD H can also be 
substantially improved. However, the gluon GPD E G is without further assumptions not 
accessible. To obtain a probabilistic interpretation, a model dependent extrapolation to the 
zero-skewness GPD has to be performed. We adopt in our studies the popular GPD model 
hypothesis that the skewness effect is t-independent, however, we also pointed out that this 
hypotheses may not hold under evolution. Certainly, this extrapolation may be considered 
as the largest theoretical uncertainty. Concerning the extrapolation uncertainties in the 
unmeasured —t region, we found that the extrapolation into the small —t region is well 
under control if the transverse resolution is of the order of 200 MeV. This allows to resolve 
the transverse distribution of partons up to 1.5 fm. In the case that the CFFs decrease 
strongly with increasing — £, e.g., exponentially, the experimentally accessible range, which 
overlaps with DVCS kinematics, is sufficient to provide an image of the sea quark GPDs H 
and E and also for the gluon GPD H G . If this will not be the case, the imaging procedure 
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may not be under control if one simply neglects the non-accessible large —t region. Under 
these circumstances one may increase experimentally the —t range together with Q 2 or 
employ mathematical extrapolation methods, which we did not explore here. 

Let us also emphasize that the revealing of the CFF £ in the small-xe region is of 
more general interest, since it is loosely related to the problem whether the "pomeron" 
coupling can flip the spin of the proton. In the partonic language it is related to the 
question of whether sea quarks and gluons carry a non- vanishing gravitomagnetic moment 
or their angular momentum is simply given by half of their momentum fraction. Certainly 
the phenomenological access to this problem suffers from the uncertainties of extrapolation 
to the forward kinematics; however, the experimental measurement of CFF £ can shed 
light on these questions, which at present can be hardly addressed with lattice gauge field 
simulations. 

In summary, the proposed high-luminosity Electron Ion Collider, combined with its 
designated detector, would be an ideal apparatus for precise measurements of exclusive 
channels in both electron-proton and electron-nuclei scattering. Besides hard exclusive 
vector meson and photon electroproduction, one might address other exclusive channels, 
too. In particular, utilizing Monte-Carlo simulations and GPD fitting routines, we have 
shown the potential of such experiments for GPD phenomenology and the ability to obtain 
from such measurements the spatial distributions of sea quarks and gluons. 
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A Updates to the MILOU code 

The MILOU code has been modified from its original version and it is currently maintained 
at BNL n . The updates to the code mainly include bug fixing together with an improved 
output. The most relevant updates are the following: 

• Bug fixed in the FORTRAN common blocks. Now they preserve the random 
seeds set in the cards and there is no need for recalculating the integral every event 
generation. 

• The correct ALLM parametrization for the F2 structure function has been 
implemented. This is relevant when running MILOU using the option for the 
Frankfurt-Freund-Strikman (FFS)[151] model, which computes the complex DVCS 
amplitude to LO and is not based on GPDs. Formerly, a wrong implementation of 
the ALLM parametrization caused a disagreement between the NLO GPD based and 
the FFS models. The correct ALLM is now taken from [152], and the agreement 
between the two models and with the predictions from the GenDVCS [153] Monte 
Carlo (also using FFS) at HERA energies is satisfactory. 

• A new output format. Beside the original output in the form of a PAW n-tuple 
[154], a new output has been implemented in the form of a Pythia-like ascii format 
text file, in the same standard as other MCs used at EIC. A detailed description of 
the new output can be found on the web-page in footnote 11. 

• Simulation of harmonics. In calculating the beam charge asymmetry, a functional 
form for the cos(cf)) harmonic was formerly hard coded. Now the code points to the 
correct values from Freund/McDermott model. 

• Simulation of the interference term, it is now properly set to the values expected 
from the Freund/McDermott model at NLO, without the twist-three contribution. 



1 A\\ information can be found at the dedicated page: https://wiki.bnl.gov/eic/index.php/MILOU 
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